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EXBCUnVB  SUMMARY 


We  are  evaluating  the  potential  impacts  of  Libby  Reservoir 
operation  on  the  fishery  in  Libby  Reservoir.  The  sampling  program 
has  been  tested  and  modified  to  provide  data  for  developing  an 
understanding  of  how  reservoir  operation  impacts  the  reservoir 
fishery.  This  sampling  program  was  reviewed  by  other  reservoir 
researchers  and  found  to  be  adequate.  Temperature  appears  to  be  an 
important  variable  influenced  by  reservoir  operation  which 
regulates  fish  and  fish  food  production  and  distribution. 

Horizontal  distributions  of  Sal  mo  spp.  showed  seasonal  move- 
ments towards  near-shore  zones  when  temperatures  allowed  and 
apparent  concentrations  in  upper-reservoir  areas  during  the  spring 
and  fall.  Sampling  fish  macro-habitat  preferences  within  the 
reservoir  proved  to  be  difficult,  but  limited  sampling  showed 
little  preference  exhibited  for  any  specific  cover  or  substrate 
types.  Juvenile  Salmo  spp.  were  found  in  both  near-shore  and 
limnetic  zones  with  smaller  (<120  mm)  juveniles  found  only  in  near- 
shore zones  and  larger  (120  to  200  mm)  juveniles  found  in  the 
limnetic  zone.  Kokanee  were  distributed  primarily  in  the  limnetic 
zone,  except  during  the  spring  and  fall  when  kokanee  were  captured 
in  near-shore  zones.  The  fall  distribution  of  kokanee  in  near- 
shore zones  was  believed  related  to  pre-spawning  movements  of 
adults  seeking  spawning  areas.  We  are  unsure  of  the  reasons  for 
spring  movements  into  near-shore  zones,  but  believe  it  may  repre- 
sent a food  seeking  response.  NOngame  species,  particularly  pea- 
mouth,  appeared  to  move  into  near-shore  zones  as  water  temperatures 
warmed  during  the  spring,  remained  in  near-shore  zones  during  the 
summer,  and  moved  out  of  near-shore  zones  in  the  fall.  Vertical 
distributions  of  Salmo  spp.  and  kokanee  appeared  to  be  regulated  by 
temperature,  but  thermal  maximums  were  occassionally  challenged, 
probably  in  order  to  feed. 

Purse  seine  population  estimates  for  Salmo  spp.  showed 
moderate  precision  for  limnetic  zones  and  poor  precision  when 
expanded  to  include  near-shore  zones.  Spring  estimates  were 
positively  biased  because  of  non-random  sampling  to  test  the  purse 
seine's  effectiveness.  An  estimated  2.5  million  kokanee  (±8%)  were 
present  during  August,  1984  based  on  hydroacoustic  sampling.  The 
presence  of  a single  (1983)  year  class  simplified  estimation. 

Kokanee  grew  rapidly  and  averaged  185  mm  after  their  first 
summer  in  the  reservoir.  Their  growth  slowed  during  the  winter  and 
averaged  only  5 mm  from  November  through  i^ril,  but  increased  again 
in  May  to  average  15  mm  per  month  through  the  summer,  ^proximate- 
ly  30  to  40  percent  of  the  male  kokanee  and  4 percent  of  the  female 
kokanee  were  determined  to  be  ripe  (based  on  August  1984  samples) 
and  probably  spawned  during  the  fall  of  1984.  These  fish  spawned 
at  age  I,  after  spending  only  two  growing  seasons  in  Libby 
Reservoir . 
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A total  of  901  adult  and  5,597  juvenile  Salmo  spp.  were  tagged 
during  1983  and  1984.  Returns  to  date  have  been  11  and  0.3  percent 
for  tagged  adults  and  juveniles,  respectively.  Movement  of  both 
juvenile  and  adult  trout  within  the  reservoir  and  between  the 
reservoir  and  its  tributaries  appeared  to  be  independent  of  tagging 
location.  A possible  down-reservoir  trend  was  noted. 

A total  of  354  (±3)  adult  cutthroat  trout  ascended  Young 
Creek  to  spawn.  The  maximum  estimated  spawning  mortality  was  25  to 
30  percent.  At  least  nine  percent  of  the  adult  cutthroat  trout 
which  spawned  during  1984  were  repeat  spawners  from  1983  based  on 
tag  returns. 

Surveys  of  tributaries  indicated  six  tributaries  contained  at 
least  one  reach  with  excellent  spawning  habitat.  One  of  these 
streams,  Barron  Creek,  is  known  to  be  unavailable  to  spawning 
adults  because  of  an  impassable  barrier  located  within  the  drawdown 
zone. 


Zooplankton,  particularly  Daohnia  spp.  and  Leotodora  spp.,  and 
terrestrial  insects  dominated  the  stomach  contents  of  Salmo  spp. 
during  the  summer  of  1983.  Kokanee  ate  Daphnia  spp.  almost 
exclusively.  These  fish  species  seemed  to  select  for  the  larger 
Dar^nia  spp.  present  in  the  reservoir.  Zooplankton  densities 
peaked  during  the  spring  with  two  lesser  peaks  evident;  one  in  the 
fall  and  one  in  mid-winter.  Daohnia  spp.  densities  usually  made  up 
20  to  30  percent  of  total  zooplankton  densities.  Length  distri- 
butions of  Daohnia  spp.  illustrated  that  length  composition  changed 
seasonally  with  more  large  Daohnia  spp.  found  during  the  fall. 
Zooplankters  were  concentrated  in  the  upper  12  m of  the  water 
column  during  the  day  with  higher  densities  of  Daohnia  spp. 
generally  found  within  the  euphotic  zone. 

Density  and  biomass  of  benthic  macroinvertebrates  were 
daviously  a function  of  how  frequently  the  reservoir  bed  was  wetted 
with  the  highest  densities  and  biomasses  observed  in  permanently 
wetted  portions. 

Abundance  of  macroinvertebrates  on  the  reservoir's  surface  was 
similar,  although  patchy,  in  near-shore  and  limnetic  zones. 

Diptera  adults  made  up  a significant  portion  of  the  total  densities 
of  macroinvertebrates  observed  during  the  spring  with  terrestrials 
making  up  a progressively  larger  percentage  of  the  total  as  summer 
progressed . 

Strategies  for  developing  predictive  models  to  describe 
relationships  between  reservoir  operation  and  the  reservoir  fishery 
will  concentrate  on  development  of  conceptual  models  first  and  then 
attempt  to  develop  quantitative  models  using  thermal  predictive 
models,  trophic  transfer  models,  and  population  simulation  models. 
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IMIRODOCnON 


Libby  Dam  was  constructed  on  the  Kootenai  (spelled  Kootenay  in 
Canada)  River  as  part  of  an  international  Columbia  River  Treaty 
between  the  United  States  and  Canada  to  provide  hydroelectric  power 
and  flood  protection  for  the  Kootenai  and  Columbia  River  basins 
(Columbia  River  Treaty  1961).  Construction  began  in  1966, 
impoundment  was  first  achieved  on  21  March  1972,  and  a full  pool 
elevation  of  2,459  feet  was  first  reached  in  July  1974. 

In  1980,  Congress  passed  the  Pacific  Northwest  Electric  Power 
Planning  and  Conservation  Act  (Public  Law  96-501)  which  created  the 
Northwest  Power  Planning  Council  (Council)  and  directed  it  to 
"promptly  develop  and  adopt.  . . a program  to  protect,  mitigate, 
and  enhance  fish  and  wildlife,  including  related  spawning  grounds 
and  habitat,  on  the  Columbia  River  and  its  tributaries."  The 
following  recommendations  by  the  Montana  Fish,  Wildlife  and  Power 
Ad  Hoc  Committee  (compiled  by  Graham  et  al.  1982)  were  adopted  by 
the  Council  as  part  of  that  program: 

1)  Except  in  years  of  extreme  runoff  (defined  as  twentieth 
percentile  or  higher  flow)  drawdown  for  power  purposes 
shall  not  exceed  90  to  110  feet  at  Libby  Reservoir 
[804(b) (1)]; 

2)  Bonneville  Power  shall  fund  research  to  develop  operating 
procedures  for  establishment  of  reservoir  levels  necessary 
to  maintain  or  enhance  fisheries  [804(b)(3)];  and 

3)  The  U.S.  Army  Corps  of  Engineers  (USACOE)  shall  develop 
operating  procedures  for  Libby  Dam  to  ensure  that  suffi- 
cient flows  are  provided  to  protect  the  resident  fish  in 
the  Kootenai  River  and  Lake  Koocanusa  (Libby  Reservoir) 
and  that  in  the  event  of  a conflict  between  maintaining 
the  minimum  flows  [804(a)(f)j  and  maintaining  reservoir 
levels  [804(b)(1)],  the  USACOE  shall  consult  with  Montana 
Department  of  Fish,  Wildlife  and  Parks  (MDFWP)  to  deter- 
mine which  requirement  shall  be  preferred  (Northwest  Power 
Planning  Council  1982). 

This  study  was  initiated  May  1983  to  determine  how  operations 
of  Libby  Dam  impact  the  Libby  Reservoir  fishery  and  suggest  ways  to 
lessen  those  impacts.  The  first  year  of  the  study  was  spent 
testing  and  selecting  methodologies  for  collecting  data  and  using 
these  data  to  design  a sampling  strategy.  We  have  finalized  our 
study  objectives  to  provide  the  following: 

1)  Quantify  available  reservoir  habitat. 

2)  Determine  the  abundance,  growth,  and  distribution  of  fish 
within  the  reservoir  and  potential  recruitment  of 
salmonids  from  Libby  Reservoir  tributaries  within  the  U.S.. 
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3)  Determine  the  abundance  and  availability  of  food  organisms 
for  fish  in  the  reservoir. 

4)  Quantify  fish  use  of  available  food  items. 

5)  Develop  relationships  between  reservoir  drawdown  and 
reservoir  habitat  for  fish  and  fish  food  organisms. 

6)  Estimate  impacts  of  reservoir  operation  on  the  reservoir 
fishery. 

Operation  of  Libby  Dam  for  flood  control  and  generation  of 
hydroelectric  power  results  in  an  annual  drawdown  and  refill  of 
Libby  Reservoir.  The  USACOE  operates  Libby  Reservoir  to  reach  a 
full  pool  in  July,  begins  drafting  the  reservoir  in  the  fall, 
reaches  a minimum  pool  elevation  in  late  March,  and  begins 
refilling  the  reservoir  in  the  spring.  These  annual  fluctuations 
of  Libby  Reservoir  water  levels  constantly  change  the  amount  and 
quality  of  reservoir  habitat  available  to  fish  and  fish  food 
organisms.  This  annual  expansion  and  reduction  of  living  space 
probably  results  in  increased  competition  for  limited  resources 
necessary  for  survival  (i.e.  food,  preferred  habitat).  Coupled 
with  this  increased  competition  is  the  potential  for  lowered  food 
production  caused  by  the  reduction  in  the  water's  surface  area 
which  receives  sunlight  energy,  the  reduced  availability  of 
nutrients  from  littoral  sources,  and  the  contraction  of  the  reser- 
voir's water  surface  away  from  productive  terrestrial  food  sources. 

We  hypothesize  that  reservoir  operation  may  affect  the 
reservoir  fishery  by: 

1)  controlling  the  amount  of  reservoir  surface  area  which 
collects  incoming  solar  energy  and  terrestrial  insects; 

2)  controlling  the  quantity  and  quality  of  habitats  available 
to  zooplankton  (volume  of  water),  benthic  invertebrates 
(wetted  reservoir  bed),  and  fish; 

3)  exposing  barriers  to  upstream  migrating  adults  located 
within  the  fluctuation  zone  of  spawning  tributaries 
resulting  in  the  loss  of  spawning  habitat  to  spring 
spawners; 

4)  flushing  j^ytoplankton,  zooplankton,  and  fish  out  of  the 
reservoir  through  Libby  Dam; 

5)  dewatering  burbot  spawning  areas  after  spawning  occurs; 

6)  eliminating  "littoral"  zone  vegetation;  and 

7)  modifying  the  thermal  regime  within  the  reservoir. 
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DESCRIPnON  OF  STOEY  AREA 


Libby  Reservoir  is  located  in  northwest  Montana  (Figure  1).  A 
detailed  description  of  the  drainage  basin  was  presented  by  Bonde 
and  Bush  (1975)  and  Woods  (1982).  We  described  the  f^ysical 
environment  and  biotic  community  in  our  first  annual  report 
(Shepard  and  McMullin  1984).  A maximum  drawdown  of  27.1  m 
(reservoir  elevation  of  2370.3  feet)  was  reached  on  15  Z^ril  1984. 
The  reservoir  reached  a full  pool  elevation  of  2459  feet  on  30  July 
1984  and  drafting  below  this  full  pool  level  began  on  22  August 
1984  resulting  in  the  reservoir  being  held  at  full  pool  for  only  23 
days  (Figure  2).  The  reservoir  was  held  at  full  pool  for  a much 
shorter  time  during  1984  than  any  previous  year  the  reservoir  was 
filled  (Table  1). 
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•BANFF 


Figure  1 . Map  of  the  Kootenai  River  Basin  showing  the  location 
of  Libby  Reservoir  (Lake  Koocanusa) (from  Woods  and 
Palter  I982). 
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RESERVOIR  ELEVATION  |ftl 


Pigu;:’e  2.  Monthly  changes  in  water  surface  elevations  of  Libby 
Reservoir  during  I982  to  1984. 
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I^ble  1.  Lake-fill  time  (yrs.),  hydraulic-residence  time  (yrs.),  maximum  drawdown,  number  of 
days  held  at  full  pool,  and  maximum  reservoir  elevation  for  Libby  Reservoir  by  year 
from  1972  through  1984. 


Lake-fill  time  (yrs)  Hydraulic-residence  (yrs)  Maximum  No.  days  Max.  pool 


Monthly 

monthlv 

drawdown 

at  full 

elevation 

Year 

annual 

mean 

min. 

max. 

annual 

mean 

max. 

min. 

ft. 

m.  ■ 

pool 

(ft.) 

1972 

0.14 

0.17 

0.04 

0.52 

0.14 

0.14 

0.02 

0.37 

1973 

0.22 

0.40 

0.10 

0.88 

0.33 

0.49 

0.11 

1.29 

230 

70.3 

1974 

0.28 

0.61 

0.09 

1.28 

0.29 

0.33 

0.13 

0.67 

153 

46.6 

1975 

0.37 

0.63 

0.11 

1.13 

0.41 

0.78 

0.10 

2.66 

172 

52.4 

0 

2454 

1976 

0.38 

0.71 

0.13 

1.54 

0.38 

0.55 

0.13 

1.56 

152 

46.2 

61 

2459 

1977 

0.64 

0.93 

0.26 

1.64 

0.50 

0.59 

0.24 

1.42 

100 

30.5 

0 

2414 

1978 

0.43 

0.75 

0.18 

1.33 

0.48 

0.63 

0.24 

1.28 

129 

39.3 

44 

2459 

1979 

0.66 

1.08 

0.22 

1.78 

0.62 

0.97 

0.22 

2.08 

95 

28.9 

0 

2451 

1980 

0.52 

0.94 

0.17 

1.47 

0.58 

0.78 

0.29 

2.07 

106 

32.3 

65 

2459 

1981 

0.33 

0.89 

0.12 

1.77 

0.41 

0.59 

0.23 

1.29 

no 

33.5 

58 

2459 

1982 

0.46 

0.84 

0.11 

1.53 

0.46 

0.49 

0.24 

0.89 

117 

35.7 

61 

2459 

1983 

0.50 

0.81 

0.13 

1.96 

0.51 

0.61 

0.22 

1.57 

111 

33.8 

54 

2459 

1984^ 

0.88 

0.16 

1.54 

0.81 

0.26 

1.53 

89 

27.1 

23 

2459 

Average 

(Post  1974)  119  52.3 


^ Through  October 
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HE^TBODS 


Metric  units  are  used  throughout  this  report  except  for 
reservoir  elevation,  reservoir  volume,  reservoir  area,  and  stream 
discharge  which  are  reported  in  feet  above  mean  sea  level,  acre- 
feet,  acres,  and  cubic  feet  per  second  (cfs),  respectively.  We  are 
using  this  convention  because  these  are  the  units  used  by  water 
managers , 


SEASONS 

We  stratified  the  year  into  four  seasons  based  on  reservoir 
operation: 

1)  Winter  (January  through  March)  - when  the  reservoir  is 
evacuated  for  flood  control  and  power  production; 

2)  Spring  (April  through  June)  - when  the  reservoir  is 
refilled; 

3)  Summer  (July  through  September)  - when  the  reservoir  is 
held  at  full  pool;  and 

4)  Fall  (October  through  December)  - when  drafting  of  the 
reservoir  begins  for  power  production. 


RESERVOIR  HABITAT 


Stratification 

The  reservoir  was  segregated  into  three  geographic  areas  based 
on  reservoir  morphometry,  effects  of  reservoir  drawdown,  and 
political  boundaries  (Figure  3).  Within  each  of  these  geographic 
areas  we  placed  permanent  sample  buoys  (the  USGS  maintains  a sample 
buoy  near  Tenmile  Creek  which  we  used)  where  water  quality  and 
zooplankton  sampling  were  conducted  (Figure  4). 

In  addition  to  these  permanent  buoys,  eight  to  12  transects 
were  established  between  recognizable  landmarks  in  each  of  the 
three  geographic  areas.  These  transects  were  subdivided  into  east, 
west,  and  mid-reservoir  stations  for  random  sampling. 

We  have  adopted  the  term  "near-shore"  zone  because  a true 
littoral  zone  does  not  develop  in  a fluctuating  reservoir  environ- 
ment. We  have  defined  the  near-shore  zone  to  be  that  portion  of 
the  reservoir  within  100  m of  the  shoreline.  That  portion  of  the 
reservoir  not  included  in  the  "near-shore"  zone  was  designated  as 
the  limnetic  zone. 

Vertical  layers  of  the  water  column  were  defined  using 
measurements  of  light  penetration,  water  temperature  (°C), 
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RESERVOIR  LATITUDINAL  PROFILES 
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Libby  Reservoir  showing 
pool  elevations. 
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4.  Map  of  Libby  Reservoir  shorn ag  sample  sites  and 
location  of  fish  traps. 
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dissolved  oxygen  (mg.liter”^) , pH,  and  conductivity  (umohs.cm"^) . 
Temperature,  dissolved  oxygen,  pH,  and  conductivity  profiles  were 
measured  with  a Martek  Mark  V digital  water  quality  analyzer  and 
light  intensity  (foot  candles)  profiles  were  measured  with  a 
Protomatic  photometer  at  permanent  buoys.  The  above  measurements 
were  made  biweekly  from  May  through  OctdDer  and  monthly  from 
November  through  April.  During  the  winter  of  1983-84  measurements 
could  not  be  made  in  the  Canada  area  from  December  through  June 
because  of  ice  and  dewatering  of  this  upper  portion  of  the 
reservoir.  The  Rexford  area  was  ice-covered  during  much  of 
December  1983. 

Depth  integrated  measurements  were  taken  at  the  surface,  at 
one  meter,  then  at  every  two  meters  down  to  15  m,  at  every  three 
meters  down  to  60  m,  and  at  every  five  meters  down  to  95  m or  the 
bottom.  Sampling  was  done  according  to  methods  used  by  the  United 
States  Geological  Survey  (USGS)  which  also  uses  Martek  Mark  V 
meters  (Greeson  et  al.  1977).  Calibrations  were  done  in  the  lab- 
oratory following  the  manufacturer's  instructions  prior  to  and 
following  field  measurements.  Accuracy  reported  by  the  manufac- 
turer were:  temperature  ±0.}Pc,  dissolved  oxygen  ±0.03ppm, 
conductivity  ±10.0  umhos.cm”^,  and  pH  +0.1  unit.  Incident  light 
was  recorded  above  the  water's  surface  and  at  one  meter  intervals 
to  a depth  of  30  m or  at  the  point  at  which  light  intensity  was  one 
percent  or  less  of  surface  light,  defined  as  the  lower  boundary  of 
the  euphotic  zone  (Greeson  et  al.  1977). 

Isopleths  of  temperature,  pH,  dissolved  oxygen,  and  conduct- 
ivity by  sample  station  across  time  and  by  date  across  sample 
stations  were  plotted  using  the  USGS  program  STAMPEDE. 


Quantity  of  Available  Reservoir  Habitat 

Contour  maps  of  the  area  impounded  by  Libby  Dam  (U.S.  Army 
Corps  of  Engineers,  Seattle  District,  File  Number  E53-1-154,  Sheets 
1-37,  1 inch  = 400  feet,  10-foot  contour  interval,  1972;  and 
British  Columbia  Ministry  of  Environment,  Map  Production,  Surveys 
and  Mapping  Branch,  Drawings  M-249-C,  Sheets  1-63,  1 inch=200  feet, 
5-foot  contour  interval,  1969)  will  be  digitized  using  a Bausch  and 
Lomb  digitizer  (Model-7048,  Huston  Instruments)  connected  to  a 
Discovery  computer.  Each  10-foot  contour  interval  will  be  entered 
by  geographic  reservoir  area.  U.S.  Army  Corps  of  Engineers 
digitized  data  for  the  above  maps  will  be  used  to  proof  our  data 
entry  and  edit  computer  data  files  of  digitized  data.  Water 
surface  area,  water  volume,  wetted  reservoir  bed  area,  and  shore- 
line length  will  be  calculated  using  a program  (GEOSCAN)  developed 
by  MDFWP  (Lonner  and  Paxton,  in  prep.). 
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Reservoir  Habitat  Happing 


During  March  1984  when  the  reservoir  elevation  was  2372  (ft) 
or  87  feet  (26.5  m)  below  full  pool,  visual  surveys  were  conducted 
covering  approximately  39.6  km  of  shoreline  in  the  Tenmile  area  and 
approximately  33.4  km  in  the  Rexford  area.  We  mapped  the  dewatered 
reservoir  bed  (from  full  pool  elevation  down  to  the  2380  foot 
contour)  and  component  maps  were  drawn  identifying  substrate  size 
classes  and  physical  cover  types  (Table  2).  Densities  of  cover 
types  were  expressed  as  the  number  of  cover  items  per  100  m of 
shoreline  (for  trees  and  logs)  or  per  100  m^  of  reservoir  bed  (for 
stumps) . 


Tributary  Stream  Surveys 

All  tributaries  to  Libby  Reservoir  within  the  United  States 
were  evaluated  for  their  potential  fish  production.  All 
tributaries  were  identified  on  USGS  quad  sheets  (scale  1:24,000) 
and  ranked  according  to  stream  order.  First  order  headwater 
segments  of  tributaries  were  assumed  to  have  little  fish  production 
potential  and  no  further  evaluations  were  conducted  on  these  small, 
steep  headwater  tributary  segments.  Tributary  streams  were  broken 
into  relatively  homogeneous  reaches  based  on  channel  gradient, 
valley  characteristics,  and  relative  flow  contribution.  This 
resulted  in  110  reaches  in  64  named  tributaries. 

Each  reach  was  surveyed.  Surveys  consisted  of  either  a 
cursory  field  exam  (if  the  reach  was  found  to  have  limited  fish 
production  potential)  or  a complete  stream  survey  in  a one 
kilometer  long  section  representative  of  the  reach  using  a 
modification  of  methods  presented  by  MDFWP  (1983)  (Appendix  A). 
Environmental  variables  sampled  included: 

a)  streambed  composition; 

b)  flood  signs; 

c)  average  depth; 

d)  mximum  depth; 

e)  channel  width; 

f)  wetted  width; 

g)  debris  accumulation  and  stability; 

h)  instreani  and  overhead  cover  (percent) ; 

^90  measurement; 

j)  compaction  and  imbeddedness  of  streambed; 

k)  composition  of  habitat  types  (pools,  riffles,  runs, 
pocketwaters,  side  channels) ; and 

l)  surface  area  of  suitable  spawning  material  in  the 
streambed. 

Measurements  or  visual  estimates  of  the  above  variables  were 
made  at  thirty  randomly  selected  transects  within  each  survey 
section.  Suitable  spawning  habitat  was  defined  as  areas  of  the 
streambed  larger  than  0.5  m^  predominated  by  material  between  2 and 
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Table  2.  Substrate  size  classes  and  j^iysical  cover  types  used  for 
mapping  the  drawdown  zone  in  Libby  Reservoir. 


Substrate  Size  Classes 

Size 

Class 

Sand  and  silt 
Gravel 
Cobble 
Boulder 
Bedrock 
Smooth 
Irregular 

<2.0  mm 
2.0  to  64  mm 
64  to  256  mm 
>256  mm 

Phvsical  Cover  Tyoes 

Abundance 
Exoressed  As 

Trees 

Stumps 

Conplex  rock  structure 

Manmade 

None 

Number  per  100  m of  shoreline 
Number  per  100  m^  of  reservoir  bed 
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64  mm  in  diameter.  The  total  area  of  suitable  spawning  habitat  in 
each  one  kilometer  section  surveyed  was  measured.  All  potential 
barriers,  sediment  sources,  and  other  unusual  landtype  features 
dDserved  within  the  survey  section  v/ere  noted. 


Enhancenent  of  Reservoir  Habitat 

One  impact  resulting  from  a fluctuating  reservoir  environment 
is  the  loss  of  terrestrial  vegetation  in  the  fluctuation  zone  and 
the  lack  of  littoral  zone  development  due  to  water  level  fluctu- 
ations (Ploskey  1982).  Potential  effects  attributed  to  the  loss  of 
near-shore  vegetation  in  Libby  Reservoir  include  the  reduction  in 
numbers  of  redside  shiners  because  of  a reduction  in  suitable 
spawning  substrates  (Huston  et  al.  1984),  increased  sedimentation 
caused  by  fluctuating  water  levels  washing  exposed  banks,  and  a 
reduction  in  the  ability  to  cycle  nutrients  through  vegetation.  We 
wanted  to  determine  the  potential  of  revegetating  the  fluctuation 
zone  in  Libby  Reservoir.  Revegetation  of  the  fluctuation  zone  was 
shown  to  be  a viable  alternative  in  other  reservoirs  (Brouha  and 
von  Geldern  1979;  personal  communication,  Del  Skeesick,  U.S.  Forest 
Service,  Portland,  Oregon) . 

Bristow  Creek  bay  and  Tobacco  River  bay  were  selected  as 
sample  sites  (Figure  4).  Unrooted  willow  (Salix  spp.)  and  red- 
osier  dogwood  (Cornus  sericea)  slips  were  planted  in  sandy  soil  in 
four  0.9  by  1.8  m plots  at  elevations  of  approximately  2454  (ft), 
2449,  2444,  and  2439  (or  depths  of  approximately  5,  10,  15,  and  20 
feet  below  full  pool)  in  April  1984.  Willow  slips  were  treated  in 
a Rootone  solution  for  approximately  two  hours  prior  to  planting. 
The  sources  for  willow  cuttings  were  an  area  near  Stone  Hill,  east 
of  the  reservoir,  and  a bottom  along  Dodge  Creek  (elevations  of 
3000  and  2750,  respectively)  for  the  Tobacco  bay  site;  and  a low 
elevation  bottom  along  Libby  Creek  (elevation  of  2360  feet)  for  the 
Bristow  bay  site.  No  dogv/ood  was  planted  in  the  Bristow  bay  site. 
Twenty-four  groups  of  two  or  three  sprigs  were  planted  in  each 
plot. 


Sedge  (Car ex  r os t rata)  sod  was  collected  from  a small  lake  bed 
near  Young  Creek.  In  May  1984,  three  pieces  of  sod  were  planted  in 
each  of  four  plots  at  both  sites,  similar  to  the  sample  design 
described  above. 

Survival  was  evaluated  in  June  1984  prior  to  innundation  by 
Libby  Reservoir.  V7e  assumed  willow  sprigs  were  alive  if  they 
exhibited  leaf  or  shoot  growth  and  sedge  survival  was  based  on  the 
presence  of  green  healthy  leaves.  Survival  will  be  assessed  again 
during  the  spring  of  1985. 
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FISH  ABUM3ANCE,  GROWTH,  At®  DISTRIBUTION 

Libby  Reservoir  presently  supports  17  species  of  fish 
(Table  3).  We  selected  westslope  cutthroat  trout,  rainbow  trout, 
hybrids  between  rainbow  and  cutthroat  trout,  kokanee,  and  burbot  as 
our  target  gamefish  species.  Northern  squawfish,  peamouth,  and 
redside  shiners  were  selected  as  our  target  nongame  fish  species. 
Information  on  distribution,  relative  abundance,  and  growth  were 
collected  for  all  fish  species  within  Libby  Reservoir,  but  our 
effort  concentrated  on  target  species  with  an  emphasis  on  gamefish 
species.  Abbreviations  used  throughout  this  report  for  the  various 
fish  species  are  as  follows:  rainbow  trout  (RB),  westslope 

cutthroat  trout  (WCT),  h^rids  betv;een  rainbow  and  cutthroat  trout 
(HB  or  WCTXRB) , members  of  the  genus  Sal  mo  spp.  (SALMO) , kokanee 
(KOK),  bull  trout  (DV),  mountain  whitefish  (MWF),  burbot  (LING), 
peamouth  (PM),  northern  squawfish  (NSO),  redside  shiners  (RSS), 
largescale  suckers  (CSU),  longnose  suckers  (FSU),  and  yellow  perch 
(YP)  . 


We  found  it  was  extremely  difficult  to  distinquish  between 
rainbow  trout,  cutthroat  trout,  and  hybrids  between  these  two 
species  using  external  morphological  characteristics.  We  based  our 
identifications  on  scale  size,  presence  or  absence  of  basibranchial 
teeth,  spotting  pattern,  and  presence  or  absence  of  the  red  "slash" 
on  each  side  of  the  jaw  along  the  dentary. 


Fish  Abundance  Trends 


Near-Shore  Zone 

Seasonal  and  annual  changes  in  fish  abundance  in  near-shore 
zones  were  assessed  using  floating  and  sinking  horizontal  gill 
nets.  These  nets  were  38.1  m long  and  1.8  m deep  and  consisted  of 
five  equal  panels  of  19,  25,  32,  38,  and  51  mm  mesh.  A floating 
gill  net  set  consisted  of  two  floating  nets  tied  end-to-end 
creating  a 76.2  m long  net.  Sinking  nets  were  set  singly.  All 
nets  were  set  perpendicular  from  shore.  Catches  were  reported  as 
the  number  of  fish  per  single  net. 


From  July  1983  to  June  1984  one  to  seven  double  floating  and 
two  sinking  gill  nets  were  set  monthly  in  each  area.  No  sampling 
was  conducted  in  the  Rexford  area  in  December  or  in  the  Canada  area 
from  December  through  June  due  to  the  presence  of  surface  ice  or 
dewatering  (Canada  area).  We  decided  to  switch  to  a seasonal, 
rather  than  monthly,  sampling  schedule  beginning  in  August  1984. 

Our  reason  for  this  change  was  to  allow  for  a more  concentrated 
effort  and  resulting  larger  sample  sizes.  Water  temperature  and 
reservoir  operation  defined  the  seasons  as  follows: 

1)  Winter  - reservoir  elevation  falling  and  reservoir 

temperature  profile  homothermous  with  surface 
water  temperatures  below  7.5°C. 
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Table  3.  Present  relative  abundance  (A=abundant,  C=common,  R=rare) 
and  abundance  trend  from  1975  to  1982  (I=increasing, 
S=stable,  D=decreasing)  of  fish  species  present  in  Libby 
Reservoir, 


Common  Name  ^ 

Scientific  name^ 

Relative 

afcaindance 

Abundance 

trend 

Gamefish  soecies 
Westslope  cutthroat 
trout 

Salmo  Clark i lewis i^ 

A 

ss/ 

Rainbow  trout 

Salmo  aairdneri 

A 

I 

Bull  trout 

Salvenlinus  confluentus 

C 

s 

Brook  trout 

Salvenlinus  fontinalis 

R 

S 

Lake  trout 

Salvenlinus  namaycush 

R 

Kokanee  salmon 

Oncorhynchus  nerka 

C 

Mountain  whitefish 

Prosooium  williamsoni 

C 

D 

Burbot 

Lota  lota 

C 

I 

Largemouth  bass 

Mic rooter us  salmo ides 

R 

D^ 

White  sturgeon 

Acioenser  transmontanus 

R 

Nonaame  ^ish  species 

Purrpk  inseed 

Leoomis  aibbosus 

R 

S / 

Yellow  perch 

Perea  flavescens 

R 

Redside  shiner 

Richardsoniusi  balteaus 

C 

d3/ 

Peamouth 

Mvlocheilus  caurinus 

A 

I 

Northern  squawfish 

Ptvchocheilus  oreaonensis  A 

S 

Largescale  sucker 

Catostomus  macrocheilus 

A 

S 

Longnose  sucker 

Catos tomus  catostomus 

C 

D 

^ From  American  Fisheries  Society  (1980) . 

^ We  adopted  the  subspecies  classification  of  Behnke  (1979) . 

^ Population  is  supplemented  with  releases  of  hatchery  origin 
fish. 

^ Kokanee  salmon  abundance  has  increased  dramatically  recently 
due  to  an  unauthorized  release  of  salmon  believed  to  originate 
from  the  Kootenay  Trout  Hatchery,  B.C. 

^ Five  white  sturgeon  were  relocated  from  below  Libby  Dam  to  the 
reservoir.  At  least  one  of  these  fish  moved  up  river  out  of 
. the  reservoir  and  two  were  reported  caught  by  anglers. 

^ Increasing  trend  for  yellow  perch  based  on  first  occurrence  in 
. recent  gill  net  catches. 

3/  Decreasing  abundance  of  redside  shiners  was  based  on  gill  net 
catches  which  capture  only  larger  (>100  mm)  individuals. 
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2) 


Spring  - reservoir  elevation  and  surface  water  temper- 
ature rising  with  surface  water  temperatures 
between  9 and  13*^C. 

3)  Summer  - reservoir  elevation  near  full  pool  and  stable 
with  surface  water  temperatures  above  17®C. 


4)  Fall  - reservoir  elevation  and  surface  water  temper- 
ature falling  with  surface  water  temperatures 
between  9 and  13°C. 

Examination  of  gill  net  catch  data  found  that  the  data 
followed  a negative  binomial  distribution.  For  comparing 
significance  of  changes  in  fish  abundance  we  can  use  the  log  (x+1) 
transformation  recommended  by  Green  (1979).  Examination  of 
probability  plotting  of  preliminary  data  adjusted  to  normal 
variates  using  the  log  (x+1)  transformation  resulted  in  a nearly 
straight  line  indicating  normality  is  a reasonable  assumption.  We 
can  detect  a difference  of  approximately  30  percent  at  an  80 
percent  level  of  significance  by  setting  approximately  20  floating 
nets  overnight  per  area  per  season.  Moyle  and  bound  (1960) 
supported  Moyle's  (1949)  contention  that  an  80  percent  probability 
is  as  high  a statistical  accuracy  that  can  be  achieved  in  practicle 
fish  population  sampling  using  passive-fishing  nets  in  lakes. 
Sinking  gill  nets  were  used  to  observe  the  vertical  distribution  of 
fish  within  near-shore  zones. 

Annual  trends  in  fish  abundance  have  been  assessed  since  1975 
by  MDFWP.  Spring  sampling  in  the  Rexford  area  with  sinking  gill 
nets  has  been  used  to  sample  the  relative  abundance  of  suckers, 
mountain  whitefish,  bull  trout,  and  burbot  (Huston  et  al.  1984). 
Fall  sampling  in  the  Rexford  and  Tenmile  (termed  Cripple  Horse  by 
Huston  et  al.  1984)  areas  with  floating  gill  nets  has  been  used  to 
sample  the  relative  abundance  of  rainbow  trout,  cutthroat  trout, 
hybrids  between  the  two  preceding  species,  redside  shiners, 
northern  squawfish,  and  peamouth  (Huston  et  al.  1984).  We 
continued  this  annual  sampling  following  criteria  described  by 
Huston  et  al.  (1984).  Spring  sampling  was  conducted  in  1975,  1976, 
1978,  1980,  1982,  and  1984.  Fall  sampling  was  conducted  in  all 
years  from  1975  through  1984  except  1977  and  1981.  Prior  to  1983 
very  few  hybrids  between  rainbow  and  cutthroat  trout  were  indent i- 
fied  although  they  were  probably  present  in  the  samples.  Investi- 
gators during  this  period  made  an  effort  to  classify  all  Sal  mo  spp. 
as  either  rainbow  or  cutthroat  trout  using  external  morphologic 
characteristics  and  the  presence  or  absence  of  basibranchial  teeth. 


Liimetic  Zone 

Vertical  gill  nets  were  used  to  assess  the  relative  abundance 
of  fish  species  in  the  limnetic  zone  of  each  area.  Replicate  over- 
night sets  of  four  45.7  m deep  by  3.7  m wide  vertical  nets  similar 
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to  those  described  by  Horak  and  Tanner  (1964)  of  mesh  sizes  19,  25, 
32,  and  38  mm  were  made  monthly  in  each  area.  These  nets  were 
marked  at  one  meter  depth  intervals.  VJe  did  not  make  replicate 
sets  in  Canada  or  in  Rexford  during  August  because  our  sampling 
buoys  were  missing.  Net  catches  were  reported  as  the  number  of 
fish  per  night  from  eight  nets.  When  only  a single  bank  of  four 
nets  were  set,  we  doubled  the  catch  for  comparative  purposes. 

A 183  m long  by  9.1  m deep  purse  seine  with  a mesh  size  of 
19  mm  was  fished  during  daylight  hours  in  the  spring  (March  through 
May)  and  fall  (October  and  November)  and  during  the  late  evening 
and  night  in  the  summer  (July).  During  the  spring  and  summer  we 
sampled  only  in  areas  where  hydroacoustic  soundings  indicated  the 
presence  of  fish  because  we  were  evaluating  the  purse  seine's 
effectiveness.  Consequently,  catches  during  these  times  were 
positively  biased. 


Vertical  Distribution 

Vertical  distribution  of  fish  was  assessed  using  vertical  gill 
nets  and  hydroacoustic  surveys.  As  vertical  nets  were  retrieved 
fish  were  removed  and  depth  of  capture  was  recorded  along  with 
species,  length,  weight,  sex,  and  state  of  maturity.  Scale  samples 
were  taken  to  determine  age. 

Hydroacoustic  sampling  was  done  across  three  permanently 
established  transects  around  the  permanent  buoys.  These  sonar 
transects  were  sampled  with  a Honda  SiTex  depth  recorder  (Model  HE- 
356A)  once  during  the  day  and  once  at  night  on  a monthly  basis  in 
conjunction  with  vertical  gill  net  sampling. 


Pqpulation  Estimates 
Purse  Seine  Estimate  of  Salmo  spp. 


The  purse  seine  described  previously  was  used  to  estimate 
numbers  of  Salmo  spp.  in  the  spring  and  fall  using  an  area-density 
estimator  (Everhart  et  al.  1975) : 


N 


A I 

a i 
i=l 


N^ 


Where: 


N = 
A = 
a = 
N^  = 


= estimated  number  of  Salmo  spp. 
total  surface  area  of  the  reservoir 
area  s£inpled  by  the  purse  seine  (0.267  ha) 
number  of  Salmo  spp.  captured  in  the  i^^ 
purse  seine  haul. 


The  estimated  variance  is: 
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a a 

a ^ N.  - ^ N. 

1=1..  . 

a(a-l) 


Assumptions  used  with  this  estimator  include  the  following: 

1)  All  the  Salmo  spp.  were  distributed  in  the  upper  6.0  m of 
the  water  column.  When  surface  water  temperatures  drop 
below  12.5^C,  vertical  nets  captured  Salmo  spp.  only  in 
the  upper  5.0  m.  If  any  Salmo  spp.  were  below  6.0  m the 
estimate  would  be  negatively  biased. 

2)  The  purse  seine  captured  all  Salmo  spp.  in  the  0.267  ha 
area  and  6.0  m depth  sampled  by  each  purse  seine  haul. 
Purse  seines  have  been  found  to  be  one  of  the  most  effi- 
cient and  least  selective  gears  for  sampling  salmonids  in 
limnetic  zones  (Hartt  1975).  If  any  Salmo  spp.  escaped 
the  purse  seine  the  estimate  would  be  negatively  biased. 

3)  The  total  surface  area  of  the  reservoir  was  known.  Area- 
capacity  tables  were  available  for  the  reservoir  (USACOE, 
1980) . 

4)  Salmo  spp.  were  distributed  uniformly  throughout  the 
limnetic  zone  of  the  reservoir,  or  if  their  distribution 
was  not  uniform,  the  sampling  distribution  was 
proportional  to  the  density  of  fish  in  different  parts  of 
the  reservoir. 

In  tl'ie  spring  39  hauls  were  completed,  primarily  in  the 
Rexford  area,  and  in  the  fall  68  hauls  were  completed  throughout 
the  reservoir.  Spring  estimates  were  positively  biased  because 
areas  sampled  with  the  purse  seine  were  selected  based  on 
hydroacoustic  locations  of  fish.  Fall  sampling  was  done  randomly. 

We  made  an  effort  to  expand  the  estimates  for  Salmo  spp.  in 
the  limnetic  zone  to  estimate  numbers  in  near-shore  zones  using  a 
correction  factor  based  on  relative  catches  in  limnetic  versus 
near-shore  floating  gill  nets.  Use  of  this  correction  factor  was 
necessary  because  densities  of  Salmo  spp.  were  higher  in  near-shore 
zones  than  in  limnetic  zones.  We  used  the  following  estimator  for 
near  shore  zones,  defined  here  as  the  portion  of  the  reservoir 
shallower  than  12.2  m (the  shallowest  depth  we  could  sample  safely 
with  the  purse  seine) : 

= (CF)Sx 

Where:  = estimate  of  Salmo  spp.  in  near-shore  zones 

corrected  using  relative  gill  net  catches; 
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CF  = correction  factor  derived  from  relative 
limnetic  versus  near-shore  floating  gill  net 
catches; 

S = number  of  sample  units  in  near-shore  zone 
(hectares  of  near-shore  zone  divided  by 
0.267  ha); 

X = mean  number  of  Salmo  spp.  captured  in  purse 
seine  hauls. 


Tlie  standard  deviation  associated  with  the  above  estimate  is: 


Hydroacoustic  Estimate  of  Kokanee 

The  number  of  kokanee  (1983  year  class)  was  estimated  using 
hydroacoustic  sampling  during  moonless  nights  in  August,  1984.  The 
kokanee  population  in  Libby  Reservoir  was  predominated  by  this  1983 
year  class  (more  than  99%).  Hanzel  (1984)  found  that  the  most 
reliable  hydroacoustic  estimates  of  kokanee  abundance  in  Flathead 
Lake,  Montana  were  made  during  moonless  nights  in  August  when 
kokanee  were  evenly  distributed  12.2  m to  36.6  m below  the  surface. 
He  attributed  this  distribution  to  the  thermal  stratification  of 
the  lake  and  the  inability  of  kokanee  to  visually  school  because  of 
darkness . 

A Honda  Sitex  depth  recorder  (Model  HE-356A)  with  a transducer 
beam  angle  of  10°  was  used.  Thirty-eight  transects  covering  76  km 
throughout  the  reservoir  were  surveyed  (Appendix  B).  All  targets 
within  each  10  m depth  interval  were  counted.  Using  vertical  gill 
net  catch  data  we  assigned  a proportion  of  the  targets  as  kokanee 
within  each  10  m depth  interval.  We  followed  methods  presented  by 
Ihorne  (1983),  but  were  unable  to  accurately  calibrate  cone  width 
in  the  field.  Vie  assumed  the  cone  width  was  consistent  with  the 
10°  beam  angle  because  we  observed  few  kokanee  above  the  10  m depth 
contour  (the  depth  zone  at  which  side  lobes  may  influence  hydro- 
acoustic target  strengths).  Since  we  were  dealing  with  a 
population  composed  almost  exclusively  of  a single  year  class  (and 
length  class)  of  kokanee,  the  problems  normally  associated  with 
variable  target  strengths  and  field  calibration  of  the  gear  were 
minimal  and  would  not  impair  our  ability  to  estimate  densities  of 
kokanee.  We  computed  the  number  of  kokanee  per  1,000  m^  for  each 
depth  interval  and  then  converted  these  number  per  volume  estimates 
to  number  per  area  estimates  using  methods  described  by  Hanzel 
(1984) . 


\/v  (5c)  = (CF) 


viler e:  V(X)  = variance  of  mean  catch  of 

Salmo  spp.  in  purse  seine  hauls 
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After  computing  the  density  of  the  1983  year  class  of  kokanee, 
an  estimate  of  the  population  was  computed  using  population 
expansion  statistics  (Snedecor  and  Cochran  1967); 


Where:  Nj^  = estimated  number  of  1983  year  class  kokanee  in 

reservoir  area  i. 

Aj^  = surface  area  of  geographic  reservoir  area  i 
(based  on  reservoir  elevation  of  2,430  feet). 
The  reservoir  was  full  (2,459  foot  elevation) 
during  sampling,  but  we  ignored  the  upper  10  m 
(30  feet) . 

Yi  = mean  density  of  1983  year  class  kokanee  in 
geographic  area  i. 

The  standard  deviation  was  estimated  by: 


H Vvar  (Yi) 

Vvar  (N^)  ^ 

\/  ai 

Where  a^^  = surface  area  sampled  in  geographic  reservoir  area 
i.  We  used  the  width  of  the  sonar  beam  at  30  m 
depth  of  5.0  m to  calculate  sampled  areas. 

Separate  estimates  were  calculated  for  each  area  and  then 
summed  across  areas  to  obtain  a total  estimate.  We  stratified  the 
reservoir  into  four  areas; 


1)  Tenmile  area  - from  the  dam  up  to  Parsnip  Creek; 

2)  Peck  Gulch  area  - from  Parsnip  Creek  up  to  the  reservoir 

bridge; 

3)  Rexford  area  - from  t±e  Reservoir  Bridge  up  to  the  U.S.- 

Canada  Border;  and 

4)  Canada  area  ~ from  the  U.S.-Canada  Border  to  the  head 

of  tiie  reservoir. 

Tiie  following  assumptions  were  made: 

1)  All  kokanee  were  distributed  below  the  10  m depth 

contour.  We  knew  a few  kokanee  were  found  less  than  10  m 
deep,  but  the  hydroacoustic  gear  discrimination  at  the  0 
to  10  m depth  interval  was  too  poor  to  identify  targets. 
The  presence  of  kokanee  above  the  10  m depth  contour 
would  negatively  bias  our  estimate. 
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2) 


We  could  apportion  hydroacoustic  targets  into  kokanee  and 
non-kokanee  using  vertical  gill  net  catches. 

3)  We  were  sampling  a 5.0  m wide  area.  Cone  width  actually 
varied  between  1.75  m at  the  10  m depth  contour  to 
10.5  m at  the  60  m depth  contour. 


Reservoir  Macro-Elabitat  Preferaice 

We  segregated  the  near-shore  zone  of  the  reservoir  using 
habitat  mapping  information  collected  on  cover  types,  substrate 
types,  and  reservoir  bed  gradient  within  the  1984  drawdown  portion 
of  the  reservoir  bed  (between  elevation  contours  2459  and  2380).  We 
returned  to  these  areas  during  the  summer  and  fall,  when  the 
reservoir  was  near  full  pool,  and  sampled  these  areas  using  boat 
mounted  electrofishing  gear  and  diver  observation.  Water  clarity 
was  so  poor  that  divers  could  not  observe  fish.  Electrofishing 
data  were  standardized  as  catch  per  unit  effort  of  electrofishing 
time. 

Relative  abundance  of  the  various  fish  species  between 
limnetic  and  near-shore  zones  were  compared  from  gill  net  data. 
Vertical  distributions  of  fish  were  compared  to  thermal  profiles 
and  vertical  distribution  of  zooplankton. 


Growth 

Fish  captured  during  all  sampling  were  measured  (total  maximum 
length)  and  weighed.  Scale  samples  were  collected  according  to 
methods  described  by  Huston  et  al.  (1984).  Ages  were  assigned 
based  on  interpreting  annulus  formation  on  scales  (Jearld  1983). 

We  estimated  growth  based  on  average  size  of  each  age  and  migration 
class  for  cutthroat  trout,  rainbow  trout,  and  hybrids  captured  in 
^ril  1984.  We  also  followed  the  growth  of  the  1983  year  class  of 
kokanee  from  November  1983  through  September  1984. 


Movement 

Fish  movement  within  the  Libby  Reservoir  drainage  was  assessed 
using  tag  return  information  and  by  trapping  several  tributaries. 
All  gamefish  captured  during  purse  seine,  electrofishing,  and 
stream  trap  sampling  were  tagged  with  either  a Floy  anchor-type  tag 
for  fish  longer  than  250  mm,  or  with  a Floy  dangler-type  tag  for 
fish  250  mm  or  smaller.  We  were  unable  to  tag  all  kokanee  because 
of  the  large  number  of  fish  handled.  We  adipose  clipped  kokanee  in 
the  spring  of  1984  and  flag  tagged  kokanee  larger  than  250  mm  with 
color  coded  tags  during  the  fall  of  1984. 

Tags  were  recovered  by  voluntary  angler  compliance  and  by  our 
fish  sampling  effort. 
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An  upstream  and  downstream  wolf-type  trap  was  operated  in 
Young  Creek  (Huston  et  al.  1984)  from  4 April  through  15  October 
1984  (Figure  4).  This  trap  was  operated  continuously  from  15  April 
through  28  July,  and  was  checked  intermittently  from  then  on. 
Downstream  box  traps  were  placed  in  Bristow  (11  June  through 
23  July),  Big  (19  June  through  25  July),  Pivemile  (18  June  through 
23  July),  Pinkham  (26  June  through  25  July),  and  Fortine  (18  July 
through  17  August)  creeks  (Figure  4).  Trap  operators  often  angled 
above  traps  to  capture  spent  adults  moving  back  down  to  the 
reservoir . 


Pood  Habits 

Stomachs  from  a representative  number  of  all  fish  species 
captured  in  gill  nets  during  1983  and  1984  were  collected  and 
emptied  into  labeled  plastic  vials  with  preservative  (Table  4). 
Laboratory  analyses  have  been  completed  for  August  1983  samples. 
Numbers  and  weights  of  each  taxonomic  group  of  food  items  in  each 
stomach  were  recorded.  Subsampling  was  used  when  analyzing 
zooplankton  and  other  abundant  small  food  items.  Wet  weights  of 
all  food  categories  but  zooplankton  were  measured  to  the  nearest 
0.01  g after  removing  excess  water  by  blotting  with  a paper  towel. 

Dry  weights  and  wet  weights  of  zooplankters  ingested  by  fish 
were  estimated  using  length-weight  regressions  (Bottrell  et  al. 
1976)  and  length-weight  tables  (Cummins  et  al.  1969).  Length- 
weight  regressions  estimated  dry  weights  (Table  5).  Dry  weights 
were  converted  to  wet  weights  using  a multiplication  factor  of  10. 
Bottrell  et  al.  (1976)  shov/ed  that  this  muliplication  factor  was 
applicable  for  wet  weights  between  10  and  300  ug,  the  range  of 
sizes  we  found  occurred  most  frequently  in  the  zooplankton  of  Libby 
Reservoir.  Cummins  et  al.  (1969)  developed  tabled  values  to 
estimate  wet  and  dry  weights  of  Leotodora  spp.  by  one  millimeter 
length  classes. 

The  majority  of  the  zooplankters  ingested  by  fish  were  frag- 
mented and  we  used  identifiable  body  parts  to  count  the  number  of 
each  genus  within  each  stomach.  We  were  forced  to  estimate  lengths 
of  Bosmina  spp.,  Diaotomus  spp.,  Eoischura  spp.,  and  Cyclops  spp. 
as  being  0.3,  0.7,  1.2,  and  0.5  mm,  respectively,  based  on  average 
lengths  found  in  the  zooplankton  population  at  the  time  food  habits 
information  was  collected.  All  dry  weights  were  converted  to  wet 
weights.  Leotodora  spp.  were  seldom  found  in  plankton  samples,  yet 
were  found  in  several  stomachs.  We  assumed  an  average  length  of 
6.0  mm  for  all  Leotodora  spp.  and  used  this  length  to  enter  length- 
weight  tables  developed  by  Cummins  et  al.  (1969)  to  estimate  the 
wet  weights  of  Leotodora.  Body  lengths  of  Daohnia  spp.  were  esti- 
mated using  measurements  of  the  post-abdominal  claw  following 
methods  presented  by  Leathe  and  Graham  (1981)  (Figure  5). 
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Table  4.  Number  of  stomachs  collected  for  food  habits  analyses  from  three  areas  of  Libby 
Reservoir  during  1983-84. 
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Table  5.  Regression  equations  used  to  estimate  dry  weights  (mg) 
of  zooplankters  from  length  data  (from  Bottrell  et  al. 
1976) . Equation:  LnW=Lna+bLnL 

where:  Ln  = natural  log; 

W = wet  weight  (mg) 

L = length  in  mm. 


Genus 


Lna  b ± 95%  C.I.  Source 


Dadinia  sod. 

1.468 

2.8292 

Bosmina  sdd. 

3.0896 

3.0395 

Cyclops  spp. 

1.9526 

2.3990 

Eoischura  sdd. 

1.9526 

2.3990 

Diaptomus  SDp. 

1.9526 

2.3990 

0.0723  Pooled  data  (Da^mia  spp. 

without  eggs,  embryos) 

0.2123  Pooled  data 

0.0854  Pooled  data  for  all 

Copepoda 

0.0854  Pooled  data  for  all 
Copepoda 

0.0854  Pooled  data  for  all 
Copepoda 
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BODY  LENGTH 


Figure  5*  Relationship  between  Daphnia  spp.  body  lengths  and 
post-abdominal  claw  lengths  measured  from  a sample 
collected  from  Libby  Reservoir  during  August,  198^. 
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An  index  of  relative  importance  (IRI)  was  calculated  to 
estimate  the  importance  of  a particular  food  item  in  the  diet 
(George  and  Hadley  1979) . The  IRI  is  the  arithmetic  mean  of  the 
number,  frequency  of  occurrence,  and  weight  of  a food  item, 
expressed  as  percentages,  in  the  diet.  IRI  values  range  from  zero 
to  100,  with  a value  of  100  indicating  exclusive  use  of  a food  item 
in  the  diet.  To  calculate  IRI's  for  insect  parts,  algae,  and 
debris  we  used  only  frequency  of  occurrence  and  weight  data. 


ABUNDANCE  AND  AVAILABILIIY  OP  PISH  FOOD  C«GANISHS 

Zoqplankton 


Standing  Crop 

Two  30.0  m (or  the  entire  water  column  when  water  depth  was 
less  than  30  m)  vertical  tows  were  made  bi-weekly  during  May 
through  October  and  monthly  from  November  through  April  in  each 
geographic  area  of  the  reservoir.  A 0.3  m diameter  Wisconsin- type 
plankton  net  was  used  for  all  tows.  One  tow  was  made  at  the 
permanent  sampling  buoy  and  one  tow  was  made  at  a randanly  selected 
point  on  an  established  transect  in  each  area.  In  June,  1984  we 
added  one  additional  random  tow  in  each  area  for  a total  of  three 
tows  per  sampling  trip  per  area.  Bottrell  et  al.  (1976)  stated 
that  variances  associated  with  replicate  samples  were  relatively 
small.  A few  replicate  tows  made  in  Libby  Reservoir  indicated  the 
same  smll  sampling  variance  existed  between  replicates,  therefore, 
no  further  replicate  sampling  was  conducted. 

A plankton  trap  similar  to  that  described  by  Schindler  (1969) 
was  used  to  sample  zooplankton  monthly  at  permanent  sample  buoys  in 
each  geographic  area.  The  trap  sampled  a 28.1  liter  volume  of 
water  and  a sample  series  consisted  of  nine  discrete  samples 
collected  at  depths  of  30,  25,  20,  15,  12,  9,  6,  and  3 m and 
immediately  below  the  water's  surface. 

All  zooplankton  samples  were  preserved  in  a solution  of  water, 
methyl  alcohol,  formalin,  and  acetic  acid.  Vertical  tow  samples 
were  diluted  in  the  laboratory  to  a concentration  in  which  five 
iTiilliliter  subsamples  contained  approximately  80  to  100  organisms. 
Schindler  trap  samples  were  concentrated  to  25  ml.  Counting  cells 
were  fabricated  out  of  Lexan  plastic  (glued  to  a glass  slide)  in 
which  a continuous  5 ml  channel  had  been  cut.  Five  individual 
5.0  ml  subsamples  were  counted  and  averaged  to  estimate  densities 
(number  per  liter)  of  zooplankters  classified  to  genus.  Carapace 
lengths  of  individual  planters  from  one  randomly  selected  5.0  ml 
subsample  were  measured.  Carapace  length  data  were  segregated  into 
0.5  mm  length  groups  for  each  genus  and  averaged  within  each  length 
group  by  genus.  Biomass  of  zooplankters  will  be  estimated  using 
the  length-weight  relationships  described  previously  in  the  "Food 
Habits"  section.  Biomass  estimates  have  not  yet  been  completed. 
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Densities  of  zooplankton  by  genus  estimated  from  vertical  tow 
samples  and  plankton  trap  sarrples  were  compared  to  evaluate  the 
efficiency  of  vertical  tows.  Schindler  (1969)  and  Bottrell  et  al. 
(1976)  considered  the  plexiglass  plankton  trap  to  be  one  of  the 
most  efficient  zooplankton  sampling  devices. 


Vertical  Distribution 

Vertical  distribution  of  zooplank ters  was  assessed  using 
plankton  trap  data  described  above. 


Benthos 

Three  replicate  samples  were  collected  seasonally  from  three 
elevational  strata  in  each  area  using  a Peterson  benthic 
dredge.  The  three  elevation  strata  were  classified  as  frequently 
dewatered  (between  full  pool  and  the  2,369  foot  contour) , occas- 
sionally  dewatered  (between  the  2,369  and  2,287  foot  contours) , and 
permanently  wetted  (below  the  2,287  foot  contour) . Benthos  sanples 
were  collected  in  October  1983,  March  1984,  and  June  1984. 

Benthos  samples  were  wet  sieved  in  the  field  using  5.6,  0.85, 
and  0.52  mm  sieves.  The  material  retained  on  the  0.52  mm  sieve  was 
preserved  and  returned  to  the  laboratory  where  all  macroinverte- 
brates were  picked  from  the  sample  and  identified  to  order  or  class 
(Diptera  and  Oligocheatea) . The  number  and  total  blotted  wet 
weight  of  each  order  (class)  were  counted  (measured)  and  densities 
were  expressed  as  number  per  square  meter  and  gran^  per  square 
meter  of  reservoir  bed  by  elevational  strata  and  geographic  area. 


Surface  Macroinvertebrates 

Macroinvertebrates  upon  the  surface  of  the  reservoir  were 
sampled  using  a net  made  of  3.17  mm  mesh  ace  bobbin  netting  tapered 
to  1.59  ram  mesh  netting  tapered  back  to  a collar.  The  mouth  of  the 
net  was  held  open  by  a rectangular  frame  measuring  1.0  m wide  by 
0.3  m high.  A removeable  plexiglass  bucket  with  a panel  of 
80  micron  Nitex  netting  was  attached  to  the  collar  at  the  cod  end 
of  the  net.  The  net  sampled  a 1.0  m wide  swath  of  the  reservoir's 
surface. 

Two  sites  were  sampled  biweekly  in  each  reservoir  area  from 
May  through  October,  and  monthly  from  November  through  April.  In 
June,  1984  we  added  one  additional  random  tow  in  each  area  for  a 
total  of  three  tows  per  sampling  trip  per  area.  Sample  sites  were 
selected  randomdy  using  established  transects  as  beginning  points 
for  the  tows.  Surface  tows  were  made  by  towing  the  net  along  the 
water's  surface  for  10  minutes  at  a fixed  rpm  which  corresponded  to 
a speed  of  1.0  m per  second.  Each  10  minute  tow  covered 
approximately  600  square  meters.  At  each  sample  site  one  tow  was 
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mde  within  100  m of  shore  (near-shore  zone)  and  one  was  made 
beyond  100  m (limnetic  zone) . 

All  macro invertebrates  were  removed  from  the  bucket  and  from 
within  the  net  at  the  end  of  each  tow  and  placed  in  a labeled  vial 
with  preservative.  All  macroinvertebrates  collected  in  each  sample 
were  identified  to  order  and  counted  at  the  laboratory.  Blotted 
wet  weights  were  measured  (in  grams)  by  order.  Densities  of 
surface  macroinvertebrates  were  expressed  as  number  per  hectare  and 
grains  per  hectare.  Paired  difference  tests  (Snedecor  and  Cochran 
1967)  were  conducted  to  determine  if  significant  differences 
existed  for  total  nun±)er  of  macroinvertebrates  per  hectare  and 
nurrtier  of  Diptera  per  hectare  between  near-shore  and  limnetic 
samples . 


RELATIONSHIP  BETWEEN  RESERVOIR  OPERATION  AND  RESERVOIR  HABITAT 

Once  we  have  established  what  types  of  reservoir  habitats  are 
important  for  producing  fish  food  organisms  and  habitats  used  by 
fish,  v/e  will  be  able  to  use  our  computer  map  data  base  with 
GEDSCAN  to  estimate  amounts  of  available  habitats  at  10  foot 
contour  intervals. 


IMPACT  OF  RESERVOIR  OPERATION  ON  FISHERY 

Our  ability  to  predict  impacts  of  reservoir  operation  on  the 
reservoir  fishery  will  depend  on  several  factors.  First,  our 
predictions  will  be  most  reliable  for  the  range  of  conditions  we 
encounter  during  the  course  of  the  study.  It  will  be  difficult  to 
reliably  predict  operational  impacts  on  the  reservoir  fishery 
unless  we  encounter  the  full  range  of  potential  operations  during 
the  study.  One  way  we  will  be  able  to  predict  operational  impacts 
will  be  through  the  use  of  predictive  models.  Secondly,  there  may 
be  a time  lag  before  reservoir  operation  effects  are  observed  in 
fishery  data.  Thirdly,  on  the  positive  side,  we  are  collecting 
information  on  all  important  trophic  variables  above  primary 
producers  and  will  be  able  to  use  these  data  to  suggest  the  best 
approach  for  predicting  operational  impacts. 


Fish  Food  Organisms 

Zooplankton  standing  crop  is  a function  of  reservoir  volume, 
primary  productivity,  thermal  regime,  flushing  rates,  and 
predation.  Reservoir  operation  dictates  reservoir  volume  and 
flushing  rates,  and  influences  primary  productivity,  thermal  regime 
and  possibly  predator  efficiency. 

Production  of  benthic  macroinvertebrates  was  related  to  the 
frequency  the  reservoir  bed  was  wetted.  The  portion  of  the 
reservoir  bed  between  full  pool  and  the  2,369  foot  elevation 
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contour  has  been  dewatered  annually  since  the  reservoir  was 
impounded  in  1972.  Ihe  portion  of  reservoir  bed  between  the  2,349 
and  2,287  foot  contours  has  been  dewatered  seven  out  of  the  last  12 
years  and  was  dewatered  last  in  1983  (Table  1) . The  reservoir  bed 
below  2,287  has  not  been  dewatered  since  full  pool  was  reached  in 
1974. 


The  density  of  macroinvertebrates  on  the  reservoir's  surface 
will  be  related  to  the  surface  area  of  the  reservoir. 


Fish 

subcontracted  with  the  USGS  to  explore  relationships 
between  reservoir  operation  and  relative  abundance,  condition,  and 
first  year  reservoir  growth  of  fish  in  Libby  reservoir  using  data 
collected  from  1974  through  1982. 

Fish  distribution  was  used  to  indicate  habitat  preference  and 
therml  preference. 


Modeling 

Four  major  types  of  models  have  been  identified  through  a 
literature  review  as  having  potential  to  predict  impacts  of 
reservoir  operation  on  the  reservoir  fishery.  A thermal  model 
developed  by  Adams  (1974)  predicts  the  thermal  regime  of  a reser- 
voir using  inflow  volumes  and  temperatures,  mixing  coefficients, 
outflow  volumes  and  depths,  and  reservoir  morphometry.  An  ad- 
vantage of  this  model  is  that  it  could  be  used  to  hold  environ- 
mental variables  constant  and  investigate  only  operational  effects. 
A fish  food  availability  versus  fish  utilization  model  (Ploskey  and 
Jenkins  1982)  has  potential  for  transferring  impacts  of  reservoir 
operation  on  fish  food  organisms  through  trophic  pathways  to 
predict  impacts  on  fish.  Habitat  suitability  models  (for  example: 
Hickman  and  Raleigh  1982)  may  be  used  to  identify  thermal  and 
chemical  preferences.  Lastly,  we  may  incorporate  a population 
simulation  model  (Serchuck  et  al.  1980)  to  follow  effects  of 
several  years  of  reservoir  operation  on  a single  target  species. 
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RESULTS 


In  this  annual  report  little  interpretation  of  results  is 
atteirpted.  Instead,  results  are  summarized  and  presented  in  graphs 
and  tables.  Detailed  summaries  are  presented  in  Appendices  C 
through  K, 


RESEEWOIR  HABITAT 


Stratification 

Surface  areas  of  each  geographic  area  within  the  reservoir  are 
impacted  differently  as  reservoir  elevation  changes  (Figure  6) . 
Storage  volumes  are  affected  similiarly.  Relative  changes  are 
largest  at  the  upper  area  (Canada)  and  smallest  at  the  lower  area 
(Tenmile) . Shoreline  lengths  and  wetted  reservoir  bed  areas  will 
be  computed  using  GEOSCAN  and  will  be  available  once  contour  data 
has  been  entered  into  a coirputer  map  data  base.  We  plan  to  have 
this  data  base  entered  into  the  corrputer  by  spring  of  1985  and 
edited  by  fall  of  1985. 

Euphotic  zone  depths  averaged  12.0,  9.5,  and  11.4  m during  the 
period  August  1983  to  October  1984  in  the  Tenmile,  Rexford,  and 
Canada  areas,  respectively  (Figure  7).  The  overall  average 
euphotic  zone  depth  for  this  period  was  11.0  m.  Typical  seasonal 
profiles  for  temperature,  pH,  and  dissolved  oxygen  (Figure  8 and 
Appendix  C)  illustrated  that  pH  and  dissolved  oxygen  levels  found 
in  Libby  Reservoir  were  generally  near  the  optimum  required  for 
cutthroat  and  rainbow  trout  (Raleigh  et  al.  1984;  Hickman  and 
Raleigh  1982) . Optimum  conductivity  levels  have  not  been  evaluated 
for  these  species. 

Vertical  layering  v/ithin  the  water  column  was  based  on  a 
combination  of  euphotic  zone  depth  and  temperature  profile  data 
because  levels  of  other  physical-chemical  parameters  were  near 
optimum  levels  reported  for  salmonids.  The  reservoir  was 
homothermous  from  October-November  through  March-April  in  the 
Tenmile  and  Rexford  areas.  The  Rexford  area  was  ice  covered  in 
Decenber  and  the  Canada  area  was  ice  covered  from  December  through 
February  when  reservoir  levels  dropped  enough  to  dewater  most  of 
that  area.  Isopleths  of  temperatures  illustrated  the  amount  of 
habitat  available  providing  optimum  tenperatures  (Figure  9) . 

Optimum  temperatures  derived  from  the  literature  for  rainbow  trout 
were  12  to  18°C  (Raleigh  et  al.  1984);  and  11  to  16*^C  for  cutthroat 
trout  (Hickman  and  Raleigh  1982) . 

The  near-shore  zone  was  defined  as  the  portion  of  the 
reservoir  within  100  m of  shore.  The  limnetic  zone  was  considered 
to  be  the  remaining  mid-water  portion  of  the  reservoir.  The  areas 
and  volumes  of  these  two  areas  will  be  estimated  using  the  GEOSCAN 
corrputer  program  and  a corrputer  map  data  base. 
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Pigiire  6.  Relationships  betvjeen  changes  in  Libby 

Reservoir  elevations  and  surface  area  of  each 
geographic  area  within  the  reservoir. 
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Figure  7.  Euphotic  zone  depths  measured  In  three  areas  of  Libby 
Reservoir  from  August ^ 19^3  through  July^  1984. 
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Pi^re  8.  Typical  seasonal  temperature,  pH,  and  dissolved  oxygen  profiles  and^ 

euphotic  zone  depths  measured  in  three  areas  of  Libby  Reservoir  during 

1984. 
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Figure  9-  Temperature  isopleths  measured  in  three  areas  of  Libby 
Reservoir  during  198^  1984.  Numbers  are  degrees 

Celsius  and  shaded  areas  represent  optimum  temperatures 
preferred  oy  cutthroat  trout  (Hickman  and  Raleigh  I982), 
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Reservoir  Habitat  Mapping 


Cover  and  substrate  types  were  mapped  (for  example, 
Figure  10) . 


Enhancement  of  Reservoir  Habitat 

All  of  the  sedge  sod  plants  survived  into  June  1984  in  both 
areas  (Table  6) . The  willow  sprigs  planted  in  the  Bristow  Bay  site 
survived  well  with  green  willow  sprigs  surviving  at  a higher  rate 
(94%  survival)  than  red  willow  sprigs  (64%  survival) . Red  willow 
and  dogwood  sprigs  planted  in  the  Tobacco  Bay  site  did  not  have  as 
high  a survival  rate  (40%  and  35%,  respectively) . 


Spawning  Habitat  in  Tributary  Streams 

Stream  reaches  were  ranked  based  on  the  amount  of  suitable 
spawning  habitat  they  contained  (Table  7) . Barron,  Bristow,  Dodge, 
Young,  Pinkham,  and  Fivemile  creeks  had  at  least  one  reach  with 
excellent  spawning  potential  (Appendix  D) . Barron  Creek  was 
unavailable  to  spawning  adult  Salmo  spp,  because  of  an  iirpassable 
barrier  located  within  the  drawdown  zone.  Big,  Cripple  Horse,  and 
Canyon  creeks  had  at  least  one  reach  with  good  spawning  potential. 
The  Tobacco  River  drainage  was  not  included  in  this  analysis 
because  habitat  survey  information  summaries  have  not  been 
conpleted. 

Kokanee  spawning  habitat  is  known  to  be  excellent  in  Kikonuin 
and  Norbury  creeks,  tributaries  located  in  Canada.  The  only  known 
spawning  habitat  used  by  significant  numbers  of  kokanee  within  the 
U.S.  is  located  in  the  Tobacco  River  from  the  reservoir  up  to  the 
town  of  Eureka.  We  assume  kokanee  have  attempted  to  spawn  along 
shoreline  areas,  but  believe  no  fry  are  produced  from  these  shore- 
line areas  because  of  dewatering  and  freezing  of  incubating 
embryos . 


FISH  ABUM5ANCE,  GROWTH,  AM)  DlSTRLBUFlOti 

Pish  Abundance 

Near-shore  Zone 

Catches  of  fish  in  floating  and  sinking  horizontal  gill  nets 
in  near-shore  zones  illustrated  the  seasonal  abundance  of  fish 
within  the  near-shore  zone  (Figures  11  and  12  and  ^pendix  E) . 

Salmo  spp.  were  captured  more  frequently  during  the  spring  in  near- 
shore net  sets  than  during  any  other  season.  This  spring  abundance 
was  probably  related  to  pre-spawning  movements  of  adults  homing  to 
spawning  tributaries  and  a food-seeking  response  to  the  increased 
availability  of  emerging  dipterans.  Salmo  spp.  were  infrequently 
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T^ble  6.  Survival  of  planted  willow,  sedge,  and  red  ozier  dogwood 
from  planting  in  May,  1984  to  June,  1984  prior  to 
inundation  by  Libby  Reservoir. 


Number  {%) 


Area 

Species 

Plot^ 

Planted 

Date 

Inventoried 

Number 

Planted 

of 

Survivors 

Bristow  Bay 

Red  Willow 

5 

04/06/84 

06/14/84 

24 

11 

(46) 

10 

It 

II 

24 

15 

(62) 

15 

II 

II 

24 

17 

(71) 

20 

n 

II 

24 

19 

(79) 

Green  Willow 

5 

If 

It 

24 

23 

(96) 

10 

n 

It 

24 

22 

(92) 

15 

II 

II 

24 

23 

(96) 

20 

It 

It 

24 

23 

(96) 

Sedge 

5 A 

05/02/84 

06/16/84 

3 

3 

(100) 

5 B 

II 

II 

3 

3 

(100) 

10  A 

It 

II 

3 

3 

(100) 

10  B 

II 

II 

3 

3 

(100) 

15  A 

II 

II 

3 

3 

(100) 

15  B 

II 

If 

3 

3 

(100) 

20  A 

II 

II 

3 

3 

(100) 

20  B 

II 

II 

3 

3 

(100) 

Tobacco  Bay 

Red  Willow 

5 A 

04/05/84 

06/23/84 

24 

16 

(67) 

5 B 

II 

II 

24 

7 

(29) 

10  A 

It 

It 

24 

8 

(33) 

10  B 

II 

II 

24 

9 

(37) 

15  A 

II 

II 

24 

10 

(42) 

15  B 

It 

It 

24 

7 

(29) 

20  A 

II 

It 

24 

8 

(33) 

20  B 

II 

II 

24 

13 

(54) 

Dogwood 

5 A 

II 

II 

24 

18 

(75) 

5 B 

II 

II 

24 

8 

(33) 

10  A 

It 

It 

24 

12 

(50) 

10  B 

II 

11 

24 

2 

(8) 

15  A 

It 

II 

24 

11 

(46) 

15  B 

It 

II 

24 

5 

(21) 

20  A 

II 

II 

24 

8 

(33) 

20  B 

II 

II 

24 

4 

(17) 

Sedge 

5 A 

05/02/84 

It 

3 

3 

(100) 

5 B 

II 

It 

3 

3 

(100) 

10  A 

II 

It 

3 

3 

(100) 

10  B 

II 

tl 

3 

3 

(100) 

15  A 

II 

II 

3 

3 

(100) 

15  B 

II 

II 

3 

3 

(100) 

20  A 

It 

II 

3 

3 

(100) 

20  B 

If 

II 

3 

3 

(100) 

^ Number  indicates  approximate  depth  at  full  pool  (ft.) 

and 

lettei 

indicates  replicate. 
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■feble  7.  Criteria  for  ranking  spawning  potential  based  on  area  of 
the  streainbed  per  kiloneter  of  stream  length  which 
contained  suitable  spawning  habitat^. 


Rank 

Square  meters  of 
suitable  spawning 
habitat^  per 
kilometer  of 
stream  lenath 

Excellent 

> 200 

Good 

> 100  < 200 

Moderate 

> 50  < 100 

Poor 

< 50 

^ For  an  explanation  of  suitable  spawning  habitat  see  the 
"Methods"  section. 
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REXFORD 


39 


Figure  11.  Monthly  catches  of  fish  (number  of  fish  per  net  night)  in  floating 
gill  nets  set  in  near-shore  zones  of  three  areas  of  Libby  Reservoi 
during  1983  and  1984. 
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caught  in  near-shore  gill  nets  set  in  the  Tenmile  and  Canada  areas 
during  the  summer,  but  surprisingly  high  nuinbers  were  captured  in 
the  Rexford  area.  We  would  expect  that  the  high  surface  water 
temperatures  observed  during  August  would  force  trout  into  deeper 
cooler  water.  The  night  of  our  August  sampling  in  the  Rexford  area 
was  extremely  windy  and  surface  temperatures  were  near  16°C 
allowing  trout  to  move  into  near-shore  zones.  Other  gamefish  most 
frequently  captured  in  near-shore  zones  included  kokanee  and  bull 
trout  which  were  caught  at  high  rates  during  spring  1984 
(Af^ndix  E)  . 

Peamouth  were  the  most  abundant  species  in  gill  net  catches 
during  the  summer.  Peamouth  appeared  to  move  to  the  bottom  of  the 
near-shore  zone  as  surface  water  temperatures  fell  during  the  fall, 
and  then  moved  out  to  the  limnetic  zone  during  winter.  They 
returned  to  near-shore  zones  in  the  early  spring  as  water 
temperatures  increased.  Redside  shiners  were  also  more  frequently 
captured  during  summer  months  in  near-shore  zones  (Appendix  E) . 
Northern  squawfish  and  largescale  suckers  mde  up  the  bulk  of 
"other  nongame  fish"  category  in  both  floating  and  sinking  gill  net 
catches . 

Annual  floating  gill  net  sampling  from  1975  through  1984 
indicated  Salmo  spp.  abundance  is  presently  lower  (1982  through 
1984)  relative  to  their  abundance  prior  to  1980  (Figure  13  and 
Appendix  F) . Abundance  of  kokanee  increased  dramatically  in  1982 
relative  to  previous  years.  The  kokanee  captured  in  1982  gill  net 
sets  were  large  (average  length  was  440  mm)  pre-spawning  (ripe) 
fish.  Floating  net  catches  of  kokanee  dropped  sharply  in  1983 
because  fish  from  the  1983  year  class  were  too  small  to  be  captured 
by  gill  nets  and  the  1981  and  1982  year  classes  were  extremely  weak 
year  classes.  Catches  of  kokanee  increased  in  1984  and  were 
comprised  primarily  of  1982  year  class  fish.  Peamouth  abundance 
has  increased  steadily  through  1982  and  appeared  to  be  leveling  off 
during  recent  years,  while  redside  shiner  abundance  dropped 
markedly  during  the  1978  to  1980  period.  Northern  squawfish 
abundance  has  fluctuated  through  the  sampling  period. 

Spring  sampling  during  the  same  time  period  with  sinking  gill 
nets  provided  no  conclusive  evidence  of  changes  in  bull  trout  or 
burbot  abundance  (Figure  13  and  Appendix  F) . The  catch  of  mountain 
v^itefish  during  the  period  1980  to  1984  was  relatively  lower  than 
that  observed  prior  to  1980.  Largescale  and  longnose  sucker 
catches  have  fluctuated  annually  with  a high  catch  of  the  former 
and  a low  catch  of  the  latter  recorded  in  1984. 


Linnetic  Zone 

Vertical  gill  net  catches  were  dominated  by  kokanee  in  the 
lower  two  areas  of  the  reservoir  (Tenmile  and  Rexford)  (Figure  14) . 
Nongame  species  dominated  the  catch  in  the  upper  (Canada)  area. 
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Figure  I3.  Annual  catches  of  fish  (number  of  fish  per  net  night)  in  floating 

gill  nets  set  during  the  fall  (top)  and  sinking  gill  nets  set  durinj 
the  spring  (bottom)  in  Libby  Reservoir  from  1975  through  198^. 
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Figure  14.  Monthly  catches  of  fish  (number  of  fish  caught  in  eight  vertical  gill 

nets  set  overnight)  in  two  areas  of  Libby  Reservoir  during  19°3  19o4. 


Catches  of  Salmo  spp.  were  highest  in  the  Rexford  area  during  the 
winter.  Catches  of  Salmo  spp.  dropped  from  May  through  July  before 
increasing  in  August.  Catches  of  bull  trout  were  consistent 
throughout  the  year.  Peamouth  dominated  the  "other”  category 
during  the  winter  and  spring.  Northern  squawfish,  largescale 
suckers,  and  mountain  whitefish  were  also  included  in  the  "other" 
category. 

Kokanee  also  dominated  the  catch  in  purse  seine  hauls 
(Table  8).  Largescale  suckers  were  the  next  most  abundant  species 
in  the  catch  followed  by  rainbow  trout  in  the  Tenmile  and  Rexford 
areas  in  the  spring  and  summer  and  the  Canada  area  in  the  fall. 
Catches  of  Salmo  spp.  were  highest  during  the  fall  in  the  Canada 
area.  Catches  of  the  Salmo  spp.  were  low,  but  second  to  kokanee, 
during  the  fall  in  the  Tenmile  area. 


Verticcil  Distribution 

Vertical  gill  net  catches  and  hydroacoustic  surveys  were  used 
to  describe  vertical  fish  distribution  and  infer  vertical  habitat 
preferences.  Vertical  gill  nets  captured  Salmo  spp.  in  waters 
cooler  than  18°C  with  cutthroat  appearing  to  prefer  slightly  deeper 
and  cooler  waters  (Figure  15  and  Z^pendix  G).  Salmo  spp.  appeared 
to  be  distributed  as  high  in  the  water  column  as  temperature 
preferences  permitted  and  generally  preferred  to  be  within  3.0  m of 
the  water's  surface  when  surface  temperatures  allowed.  Rainbow 
trout  were  occassionally  found  deeper  than  other  Salmo  spp. 

Kokanee  were  distributed  throughout  the  upper  45  m of  the 
water  column,  but  appeared  to  prefer  temperatures  cooler  than  17 ®C 
(Appendix  G).  Peamouth  were  captured  in  deep  water  (generally 
deeper  than  15  m)  throughout  the  winter  and  moved  up  towards  the 
surface  in  the  spring  as  water  temperatures  warmed.  Peamouth  were 
found  near  the  surface  during  the  summer  once  water  temperatures 
rose  above  10°C.  Bull  trout  seemed  to  select  areas  cooler  than 
14°C  and  exhibited  an  apparent  preference  for  water  10  to  12°C. 

Hydroacoustic  surveys  showed  that  targets  (fish)  were  sparsely 
distributed  across  sample  transects  during  daylight  and  were 
densely  distributed  across  transects  during  hours  of  darkness 
(Figure  16).  Vertical  distribution  from  the  surface  down  to  five 
meters  was  difficult  to  decipher  from  hydroacoustic  charts  because 
of  near  surface  interference.  Fish  distribution  below  5 meters  on 
hydroacoustic  charts  generally  corresponded  to  catches  in  vertical 
gill  nets  set  in  conjunction  with  hydroacoustic  surveys. 
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Figure  I5.  Typical  seasonal  vertical  distributions  of  fish  and  zooplankton  compared 
to  temperature  profile  and  euphotic  zone  depth  in  the  Rexford  area  of 
Libby  Reservoir  during  4 July  1984. 


Figure  16.  Hydroacoustic  charts  showing  distribution  of  fish  identified 
as  targets  during  the  day  (15OO  hours;  top  chart)  and  night 
(2230  hours;  bottom  chart)  on  2 October  1984  in  the  Tenmile 
area  of  Libby  Reservoir. 
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Population  Estimates 


Purse  Seine  Estimate  of  Salmo  spp. 

Estimates  of  Salmo  spp.  populations  within  the  limnetic  zone 
were  made  using  purse  seine  catches  during  the  spring  and  fall 
(Table  9).  Spring  estimates  were  known  to  be  high  because  we 
positively  biased  the  results  in  actively  seeking  fish  using  hydro- 
acoustic sampling  prior  to  purse  seining  to  evaluate  the  effective- 
ness of  the  gear.  We  have  no  way  of  estimating  the  magnitude  of 
our  introduced  bias. 

Fall  purse  seine  estimates  were  computed  using  two  statistical 
designs,  stratification  by  area  and  no  stratification.  Ninety-five 
percent  confidence  intervals  around  the  two  limnetic  estimates 
overlapped,  but  the  stratified  point  estimate  was  lower  and  had 
poorer  precision.  Fall  estimates  may  be  low.  However,  since 
spring  estimates  were  known  to  be  high  and  fish  remaining  in  the 
reservoir  in  the  fall  survived  through  the  summer  fishery  our  fall 
estimates  may  be  reasonable. 

Near-shore  gill  nets  captured  4.4  and  2.4  times  as  many  Salmo 
spp.  as  limnetic  sets  in  the  spring  and  fall,  respectively.  These 
values  were  used  as  correction  factors  to  expand  limnetic  purse 
seine  catches  for  estimating  numbers  of  Salmo  spp.  in  near-shore 
zones.  Using  limnetic  estimates  to  estimate  Salmo  spp.  populations 
within  near-shore  zones  resulted  in  poor  levels  of  precision 
(Table  9) . 


Hydroacoustic  Estimate  of  Kokanee 

Hydroacoustic  sampling  data  collected  during  August  1984  were 
statistically  analyzed  to  estimate  numbers  of  kokanee  (1983  year 
class).  We  estimated  there  were  slightly  more  than  2.5  million 
kokanee  (+  8%)  in  the  reservoir  at  that  time  (Table  10).  Using 
average  weight  data  the  above  number  represents  approximately 
400,000  kg  of  kokanee  or  24.5  kg  per  hectare  (21.8  lbs  per  acre). 
The  highest  densities  of  kokanee  were  found  in  the  Tenmile  and 
Rexford  areas,  while  the  Peck  Gulch  and  Canada  areas  had  signifi- 
cantly lower  densities  than  either  of  the  above  two  areas  and 
Canada  had  significantly  lower  densities  than  Peck  Gulch 
(Table  11) . 


Reservoir  Macro-Habitat  Preferences 

Vertical  and  thermal  habitats  preferred  by  fish  were  described 
previously  in  the  "Vertical  Distribution"  section.  We  attempted  to 
describe  our  present  understanding  of  areal  fish  distribution 
during  the  year  by  using  relative  horizontal  and  vertical  gill  net, 
electrofishing,  and  purse  seine  catches  of  fish  to  infer  where  fish 
preferred  to  be  during  the  course  of  a year.  When  using  passive 
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T^ble  9.  Estimated  nuirbers  of  Salmo  spp.  derived  from  purse  seine 
catches  in  Libby  Reservoir  during  spring  and  fall,  1984. 


Season 

Geographic 

Area 

Zone 

Surface 
area (acres) 

Estimate 

95%  CI^  (%) 

Spring 

Total^ 

Total 

Limnetic 

Near-shore 

Total 

21,175 

5.521 

26,696 

192,500^ 
220. 836^/ 
413,336^ 

+ 39,949  (21) 
+424.181 f 192) 
±464,130(112) 

Fall 

Total 

Total 

Limnetic 

Near-shore 

Total 

34,090 

9.522 

43,612 

42,536 

28.392 

70,928 

+ 16,491  (39) 
+ 99.082(349) 
±115,573(163) 

Tenmile 

Rexford 

Canada 

Limnetic 

Limnetic 

Limnetic 

15,875 
10,280 
8,300^  . 

6,804 

5,717 

23.795 

+ 13,336(196) 
+ 4,378  (76) 
+ 6.661  (28) 

Total 

Limnetic 

34,455^ 

36,316 

24,375  (67) 

^ 95%  Cl  is  the  nine-five  percent  confidence  interval. 

^ Total  indicates  estimate  done  for  entire  reservoir 
(No  stratification) . 

^ Ihese  estimates  were  high  because  of  positive  bias  introduced 
. by  sampling. 

^ Total  areas  do  not  match  because  we  had  to  interpolate 
between  10  foot  contours  on  stratified  estimate. 
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Table  10.  Estimated  nuirbers  of  1983  year-class  kokanee  in  Libby 
Reservoir  expanded  from  hydroacoustic  sampling  during 
August,  1984. 


Geographic 

Total  Surface  Area^ 

Area  Sampled 

Estimated 

95%  confidence 

Area 

ha 

acres 

ha 

acres 

Number 

Interval 

(%) 

Tenmile 

4,743.6 

11,721.4 

10.8 

26.8 

1,228,637 

+ 69,254 

(6) 

Peck  Gulch 

1,959.7 

4,842.3 

7.2 

17.9 

257,804 

+ 25,040 

(10) 

Rexford 

4,466.3 

11,036.1 

10.8 

26.7 

1,059,356 

+113,540 

(11) 

Canada 

-^321^ 

13.150.2 

_9^ 

22,2 

28.536 

+ 6.397 

1221 

TOTAL 

16,591.5 

40,750.0 

37.8 

93.6 

2,574,333 

±214,231 

(8) 

Total  surface  area  based  on  a reservoir  elevation  of  740.7  m 
(2,430  ft.)  (see  Methods  for  explanation). 
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Table  11.  Density  of  1983  year-class  kokanee  in  four  areas  of 

Libby  Reservoir  in  August,  1984  determined  from  hydro- 
acoustic sampling  and  analysis  of  variance  results  and 
IXincan's  new  multiple  range  test  results  (OH  1977) 
illustrating  where  densities  differ  significantly. 


■Density 


Area 
Tenmile 
Rexford 
Peck  Gulch 
Canada 


Number/acre 

104.8 

96.0 

53.2 

2.2 


MjmberAia 

259.0 

237.2 

131.4 

5.4 


Duncan's 

F-test^  test^ 

144.6**^  I 


^ Results  from  analysis  of  variance  test  ccanparing  densities  from 
all  four  areas. 

^ Duncan's  new  multiple  range  test  at  a = 0.05  level.  Solid  line 
between  connecting  areas  indicates  no  significant  difference  in 
kokanee  density  between  areas. 

^ Significant  at  p<0.001  level. 
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sampling  techniques  such  as  gill  nets,  one  must  be  careful  about 
comparing  seasonal  abundance  because  of  variable  activity  rates 
associated  with  changing  water  temperatures. 

Salmo  spp,  appeared  to  prefer  near-shore  zones  as  long  as 
temperatures  in  these  zones  were  within  tolerable  limits.  Salmo 
spp.  catches  in  near-shore  nets  were  highest  in  spring  corres- 
ponding to  the  lowest  catch  in  limnetic  vertical  nets  suggesting 
they  moved  into  near-shore  zones  during  the  spring.  Reasons  for 
this  movement  to  near-shore  zones  were  believed  to  be  related  to 
increased  food  availability  in  the  form  of  emerging  Diptera,  a pre- 
spawning movement  to  begin  searching  for  natal  spawning 
tributaries,  and  increased  metabolic  activity  in  response  to 
warming  water  temperatures.  Water  temperatures  in  near-shore  zones 
were  9 to  11®C  in  the  Rexford  area  and  6 to  11°C  in  the  Tenmile 
area  during  May  and  June.  During  the  summer,  as  water  temperatures 
in  near-shore  zones  warm,  Salmo  spp.  moved  offshore  into  deeper 
cooler  water,  but  appeared  to  occassionally  challenge  their  thermal 
habitat  preference  as  indicated  by  our  relatively  high  floating 
gill  net  catch  in  the  Rexford  area  during  August.  The  surface 
temperature  in  this  area  at  that  time  was  18.6*^C.  Salmo  spp.  moved 
into  near-shore  zones  during  the  fall  and  v/ere  found  in  both  the 
limnetic  and  near-shore  zones  during  the  winter.  We  did  see  an 
increase  in  limnetic  catches  in  vertical  gill  nets  in  the  Rexford 
area  which  may  indicate  a movement  to  the  upper  end  of  the 
reservoir  during  the  winter. 

Electrofishing  of  near-shore  zones  containing  various  habitat 
types  revealed  no  distinct  trends  in  habitat  preference  exhibited 
by  distributions  of  Salmo  spp.  (Table  12).  We  captured  more  large 
trout  (>300  mm)  and  mountain  whitefish  during  fall  sampling  than 
during  spring  sampling. 

Electrofishing  and  purse  seine  catches  of  fish  were  used  to 
describe  the  distribution  of  juvenile  Salmo  spp.  in  Libby  Reservoir 
since  juvenile  trout  were  seldom  captured  in  gill  nets.  Juvenile 
Salmo  spp.  were  captured  in  both  near-shore  and  limnetic  zones; 
however,  smaller  juveniles  (<120  mm)  were  captured  exclusively  in 
near-shore  zones  while  larger  juveniles  (120  to  200  mm)  were 
captured  exclusively  in  limnetic  zones.  We  speculate  that 
territorial  stream-rearing  juveniles  may  visually  key  on  physical 
habitat  when  they  first  enter  the  reservoir.  This  behavioral 
response  to  physical  habitat  forces  these  juveniles  to  inhabit 
near-shore  zones.  Later,  after  these  juveniles  "learn"  to  seek 
the  previously  unfamiliar  zooplankton  food  resources,  they  may  move 
offshore  into  limnetic  zones  where  food  is  abundant  and  predation 
is  reduced. 

Kokanee  preferred  limnetic  zones  throughout  the  year.  The 
only  seasons  when  kokanee  were  found  in  near-shore  zones  were  in 
the  early  spring  (/^ril)  and  the  fall.  Spring  use  of  near-shore 
zones  may  be  related  to  increased  availability  of  Diptera  coupled 
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Electrofished  during  the  c3ay 


with  lower  densities  of  Da^^ia  spp.  Movement  of  kokanee  into 
near-shore  zones  during  the  fall  may  have  been  related  to  a pre- 
spawning search  for  suitable  spawning  habitat. 

Bull  trout  used  limnetic  and  near-shore  zones,  but  were  seldom 
found  near  the  surface.  Peamouth  and  northern  squawfish  moved  into 
near-shore  zones  during  the  spring.  They  remained  near  the  surface 
in  near-shore  and  limnetic  zones  during  the  summer,  moved  to  deeper 
near-shore  habitats  during  the  fall,  and  out  into  the  limnetic  zone 
during  the  late  fall  and  winter.  Redside  shiners  were  relatively 
abundant  in  near-shore  zones  during  the  spring  and  summer,  but  were 
seldom  captured  in  either  near-shore  or  limnetic  gill  nets  during 
other  seasons. 


Growth 

Growth  of  rainbow  and  cutthroat  trout  was  evaluated  using 
April  1984  samples  of  fish  (Table  13).  Empirical  lengths  at  this 
time  should  closely  represent  annual  growth  since  this  is  approx- 
imately when  annulus  formation  occurs.  Rainbow  trout  in  the  sample 
were  predominated  by  fish  which  spent  one  (51%)  and  two  (38%)  years 
rearing  in  natal  tributaries  (migration  class  and  X2, 
respectively).  Second  year  reservoir  growth  for  these  two 
migration  classes  was  78  and  66  mm,  respectively,  based  on 
comparing  two  consecutive  year  classes.  Rainbow  from  both 
migration  classes  reached  approximately  420  mm  by  age  V.  Wild 
cutthroat  trout  were  predominated  by  migration  class  two  (66%)  and 
three  (24%)  fish.  Hatchery  origin  fish  made  up  37  percent  of  the 
total  cutthroat  trout  sample.  Second  year  reservoir  growth  was  81 
and  64  mm  for  these  two  migration  classes,  respectively,  based  on 
comparing  two  consecutive  year  classes.  Cutthroat  from  both 
migration  classes  reached  approximately  380  mm  by  age  V and  405  mm 
by  age  VI. 

Monthly  growth  of  the  1983  year-class  of  kokanee  was  followed 
from  November  1983  through  August  1984  (Figure  17).  Kokanee 
average  185  mm  in  November,  although  vertical  nets  probably 
selected  for  larger  individuals.  Growth  from  November  through  the 
following  ^ril  was  slow  and  averaged  only  5 mm.  Weight  increased 
by  an  average  of  15  grams  during  that  same  time  period.  Beginning 
in  May,  growth  increased  rapidly  and  length  increased  by  an  average 
of  15  mm  per  month  through  September,  while  weights  increased  by 
19  grams  per  month.  Using  August  sampling  we  estimated  that 
approximately  30  to  40  percent  of  the  male  kokanee  and  5 percent  of 
the  female  kokanee  spawned  during  1984  based  on  size  of  the  gonads. 
These  kokanee  were  age  I and  had  spent  only  two  growing  seasons  in 
Libby  Reservoir.  Growth  then  appeared  to  slow  again  and  the  loss 
of  some  larger  members  of  the  1983  year  class  to  the  1984  spawning 
resulted  in  a smaller  average  size  in  the  remaining  population. 
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T^ble  13.  Average  lengths  of  rainbow  and  cutthroat  trout  in  April  1984 
by  age  and  migration  class  based  on  interpretation  of  scale 
samples  taken  from  fish  captured  in  Libby  Reservoir  (Number 
of  fish  in  sample  in  parentheses). 


Species  Aae 


Migration  class 

n(%) 

II 

III 

IV 

V 

VI 

Rainbow  Trout 
y 3/ 

^2 

^3 

Average 

8 (9) 
45(51) 
33(38) 
2 (2) 

332  (7) 
320(22) 

430  (1) 
398(15) 
323(11) 

420  (4) 
389(13) 
365  (2) 

421  (4) 
419  (9) 

— 

323(29) 

369(27) 

393(19) 

420(13) 

— 

Cutthroat  Trout 

30 

4 (8) 
1 (2) 
33(66) 
12(24) 

306(13) 
308  (1) 
283  (1) 

372(10) 
356  (3) 

301(20) 

392  (7) 

382(10) 

315 

385  (2) 
379  (4) 

408  (1) 
406  (1) 

Average 

305(15) 

328(33) 

365(24) 

381  (6) 

407  (2) 

Subscripted  X indicates  age  at  which  fish  entered  Libby  Reservoir 
from  tributaries.  An  "H"  indicates  a hatchery  origin  fish. 
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Figure  I7.  Monthly  groiArth  of  the  198^  year  class  of  kokanee  in 
Libby  Reservoir  from  November,  I983  through  August, 
1984.  Vertical  bars  represent  95^  confidence  limits. 
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Young  Creek 

An  estimated  354  adult  westslope  cutthroat  trout  ascended 
Young  Creek  to  spawn  during  1984  (Appendix  H).  The  sex  ratio  of 
this  spawning  run  was  1:2.4  Spawning  mortality  was 

estimated  to  be  between  25  to  30  percent,  although  this  must  be 
considered  a maximum  estimate  because  we  had  a mink  predation 
problem  during  our  downstream  trapping.  We  tagged  a total  of  227 
adult  cutthroat  trout  moving  downstream  through  our  trap.  We 
observed  that  at  least  nine  percent  of  the  adults  which  spawned 
during  1984  were  previous  year  spawners  based  on  tag  returns.  This 
is  a minimum  estimate  because  an  unknown  amount  of  tag  loss 
occurred.  We  captured  one  adult  cutthroat  (a  ripe  male  287  mm 
long)  which  we  had  tagged  the  previous  year  as  a juvenile 
emigrating  downstream  (160  mm  long).  This  fish  grew  127  mm. 

We  trapped  a total  of  1,330  juvenile  trout  emigrating 
downstream  in  Young  Creek  (/^pendix  H).  The  majority  of  these 
juveniles  were  tagged  (982)  and  emigrated  from  mid-May  through  mid- 
July.  Two  hundred  fifty-six  of  these,  particularly  young-of-the- 
year  (57%),  emigrated  later,  primarily  during  October  and  November. 


Bristow  Creek 

A total  of  36  adult  Sal  mo  spp.  were  captured  and  tagged 
leaving  Bristow  Creek  after  spawning.  The  majority  were  cutthroat 
(55%)  and  hybrids  (33%).  Most  of  the  adults  left  by  the  end  of 
June  (/^pendix  H).  Approximately  40  percent  of  the  adult  cutthroat 
were  determined  to  be  of  hatchery  origin  based  on  eroded  dorsal 
fins.  Sex  ratios  were  1:3  (c/  :$)  for  all  Salmo  spp.  combined. 

A total  of  175  juveniles  were  trapped  and  tagged  emigrating 
from  Bristow  Creek.  The  majority  were  identified  as  cutthroat 
(94%)  and  juvenile  emigration  appeared  to  end  after  10  July, 
although  the  trap  was  operated  through  23  July  (^pendix  H) . 

Big  Creek 

A total  of  51  cutthroat  trout  (10  percent  identified  as 
hatchery  origin),  4 rainbow  trout,  and  22  hybrid  adults  were 
captured  moving  out  of  Big  Creek  after  spawning  (^pendix  H). 
Sixty-three  of  these  fish  were  tagged.  Most  of  these  fish  moved 
through  the  trap  between  4 July  and  20  July. 

A total  of  1,549  cutthroat  trout,  24  rainbow  trout,  and  278 
hybrid  juveniles  were  trapped  emigrating  downstream  in  Big  Creek 
(i^)pendix  H)  . A total  of  1,730  were  tagged.  Movement  occurred 
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throughout  the  trapping  period,  but  was  slowing  when  traps  were 
removed  on  25  July. 


Fivemile  Creek 

Rainbow  trout  adults  predominated  the  downstream  catch  of 
adult  trout  in  Fivemile  Creek  (28  fish),  followed  by  cutthroat  (21) 
and  hybrids  (19).  Downstream  movement  of  adults  was  spread  out 
over  the  trapping  period  (Appendix  H).  This  delayed  movement  may 
have  been  caused  by  fish  holding  above  our  trap  because  when 
streamflows  dropped  our  trap  was  in  a slow  run.  Approximately  38 
percent  of  the  adult  cutthroat  trout  captured  were  hatchery  origin 
fish,  similar  to  Bristow  Creek.  Sex  ratios  for  all  three  species 
were  skewed  more  towards  males  than  we  found  in  other  tributaries 
(WCT-1:0.6;  RB-1:0.8;  HB-1:1). 

The  majority  of  downstream  migrating  juvenile  trout  were 
identified  as  cutthroat  trout  (258)  with  very  few  rainbow  trout  (1) 
or  hybrids  (8)  trapped.  Juveniles  moved  primarily  in  late  June  and 
early  July  (^pendix  H).  The  absence  of  juvenile  rainbow  trout  in 
our  trap  catch  may  reflect  a problem  in  identification  or  may 
indicate  juvenile  rainbow  trout  emigrate  at  smaller  sizes  (earlier 
ages)  and  were  not  vulnerable  to  our  trap. 


Pinkham  Creek 


The  lower  portion  of  Pinkham  Creek  channel  has  a high  gradient 
and  much  of  the  streamflow  flows  subsurface.  The  downstream  trap 
in  Pinkham  Creek  was  placed  in  a poor  location  because  once  stream- 
flows  receded  the  trap  was  left  in  a slackwater  pool.  Conse- 
quently, most  of  the  adult  fish  held  above  the  trap  and  had  to  be 
captured  using  angling.  Therefore,  we  were  unable  to  use  catches 
of  fish  in  the  downstream  trap  to  describe  the  timing  of  downstream 
movement.  We  captured  and  tagged  14  cutthroat  trout,  12  hybrids,  and  five 
rainbow  trout  adults.  A total  of  55  juvenile  trout  were  captured. 

The  majority  of  these  were  identified  as  cutthroat  trout  (91%). 

Forty-seven  juveniles  were  tagged. 


For tine  Creek 

Only  three  adult  and  85  juvenile  Salmo  spp.  were  captured  in 
the  downstream  trap  placed  in  Fortine  creek  (i^pendix  H).  We  were 
unable  to  install  the  trap  until  18  July  because  of  the  size  of  the 
stream  and  volumes  of  water  discharged.  Juvenile  movement  through 
the  trap  ended  by  1 August.  We  also  captured  mountain  whitefish, 
largescale  suckers,  redside  shiners,  northern  squawfish,  and  eastern 
brook  trout  in  the  trap.  Water  temperatures  in  Fortine  Creek  were 
generally  higher  than  other  tributaries  trapped  during  1984. 
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Tag  Returns 


Movement  within  Libby  Reservoir  and  between  the  reservoir  and 
its  tributaries  was  assessed  using  voluntary  angler  tag  returns  and 
tagged  fish  recovered  by  our  sampling  gear  (^^pendix  H). 
^proximately  11  percent  of  the  901  tagged  adult  Salmo  spp.  were 
recaptured.  The  majority  (95%)  were  returned  by  anglers.  It 
appeared  adults  moved  throughout  the  reservoir  and  tributary  origin 
did  not  play  a major  role  in  reservoir  distribution  (Table  14). 
Ihere  may  be  a slight  down-reservoir  trend,  but  several  fish  moved 
up-reservoir . 

Returns  of  the  5,597  dangler  tags  tags  placed  on  juvenile 
Salmo  spp.  were  extremely  low  (less  than  one  percent).  This  low 
return  could  be  related  to  one  or  more  of  the  following:  small 

size  of  the  tags  making  them  difficult  for  anglers  to  notice;  the 
small  size  of  juveniles  preventing  them  from  making  up  a signifi- 
cant portion  of  angler's  catch;  or  tag  loss.  Using  this  limited 
information  it  appeared  juveniles  also  moved  throughout  the 
reservoir  (Table  14).  Several  were  captured  by  anglers  from  the 
tributary  in  which  they  were  tagged  shortly  after  being  tagged.  It 
was  interesting  to  note  that  both  an  adult  and  a juvenile  were 
captured  in  the  Kootenai  River  below  Libby  Dam. 


Reservoir  Food  Habits 

Food  habits  data  have  been  summarized  for  gamefish  collected 
during  the  summer  of  1983  (i^pendix  I).  IRI  values  showed  that 
zooplankton,  particularly  Daohnia  spp.  and  Leotodora  spp.,  and 
terrestrial  insects  made  up  the  bulk  of  the  diet  of  Salmo  spp. 
during  the  summer  (Figure  18).  Smaller  (<.330  mm)  Salmo  spp. 
appeared  to  prefer  zooplankton  and  Diptera  adults,  while  larger 
(>330  mm)  Salmo  spp.  seemed  to  exhibit  a preference  for  Dadinia 
spp.  and  terrestrial  insects.  Mountain  whitefish  ate  zooplankton 
comprising  all  genera  and  Diptera  adults  and  larvae.  Kokanee  ate 
Daphnia  spp.  almost  exclusively. 

Gamefish,  particularly  trout  and  kokanee,  seemed  to  select 
larger  Daphnia  spp.  as  indicated  by  comparing  Daphnia  spp.  size 
compositions  and  average  sizes  in  the  diet  with  those  sampled  from 
the  zooplankton  population  (Table  15). 


ABUNDANCE  AN)  AVAILABILZIY  OP  FISH  FOCD  0«GANISHS 

Zooplankton 


Standing  Crop 

Our  data  indicated  the  0.3  m plankton  net  was  more  efficient 
than  the  plankton  trap  (Table  16).  This  result  conflicts  with 
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Number  of  fish  taggeci. 


Figure  18.  Indices  of  relative  importance  (iRl)  for  gamefish  collected 
in  Libby  Reservoir  (all  three  areas  combined)  during 
August,  1983* 
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T^ble  15.  Size  composition  (%),  mean  lengths,  and  mean  estimated  weights  of  . Daphnia  spp. 
ingested  by  gamefish  and  available  in  the  zooplankton  population^  in  Libby 
Reservoir  during  August, 1983. 
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Available  zooplankton  as  sanpled  by  Wisconsin  tows  during  August. 
Fish  abbreviations  were  explained  in  the  "Methods"  section. 
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results  documented  in  the  literature  (Schindler  1969,  Bottrell  et 
al.  1976)  and  we  cannot  explain  why  it  occurred.  Regardless,  it 
seemed  reasonable  for  us  to  assume  plankton  net  hauls  represented 
zooplankton  densities  most  accurately,  and  densities  reported  here 
were  based  on  net  tow  data. 

Zooplankton  densities  peaked  during  July  in  the  Tenmile  area 
(14.1  per  liter).  May  in  the  Rexford  area  (18.06  per  liter),  and  in 
November  in  the  Canada  area  (27.34  per  liter)  (Figure  19).  Two 
lesser  peaks  occurred  in  both  the  Tenmile  and  Rexford  areas,  one 
higher  peak  during  S^tember  and  another  lower  one  in  January  or 
February.  Daphnia  spp.  generally  comprised  20  to  30  percent  of  the 
numbers,  while  Cyclops  sj^.  generally  made  up  the  largest 
percentage  (Appendix  J)..  Bosmina  spp.  were  most  abundant  during 
the  summer  and  present  in  low  numbers  during  the  winter. 

Length  distributions  of  Daphnia  spp.  showed  that  larger 
Daphnia  spp  (>1.5  mm)  were  most  abundant  during  the  fall  and  absent 
during  much  of  the  winter  (Figure  20).  Areal  differences  existed 
in  length  frequency  distributions  between  the  Tenmile  and  Rexford 
areas  with  few  larger  Da|dinia  spp.  seen  in  the  Rexford  area  during 
the  spring  and  early  summer.  Daphnia  spp.  smaller  than  0.5  mm  were 
not  seen  in  plankton  samples  during  the  late  fall. 


Vertical  Distribution 

Zooplankton  densities  were  higher  during  the  day  above  the 
12  m depth  contour  than  below  it  (i^pendix  G).  Densities  of 
zooplankers,  particularly  Daphnia  spp.,  appeared  to  highest  within 
the  eu^^otic  zone,  and  vertical  kokanee  distribution  appeared  to  be 
related  to  vertical  Daphnia  spp.  distribution  (Figure  15  and 
impend  ix  G) . 


Benthos 

Densities  and  biomass  of  Diptera  larvae  and  Oligocheates  were 
lowest  in  the  frequently  dewatered  portions  of  the  reservoir  bed  and 
highest  in  the  permanently  wetted  portions  in  the  Tenmile  and 
Rexford  areas  (Table  17).  The  Rexford  area  had  a higher  density 
and  biomass  of  Dipterans  than  the  Tenmile  area  when  comparisons 
were  made  within  drawdown  frequency  strata. 


Surface  Hacroinvertebrates 

Comparison  of  near-shore  versus  limnetic  surface  tows  found 
that  although  total  mean  densities  of  macroinvertebrates  were  often 
higher  on  the  limnetic  water  surface,  there  was  no  significant 
difference  between  near-shore  and  limnetic  catches  (Table  18).  The 
same  result  was  found  for  the  adult  Dipteran  component  of  the 
surface  catch.  Distribution  of  surface  macroinvertebrates  was 
found  to  be  patchy,  both  spatially  and  temporally  (Figure  21),  but 
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rCNMU 


15- 


I I CTClO^t 


Figure  19.  Densities  (number  per  liter)  of  the  four 

most  abundant  genera  of  zooplankton  in  three 
areas  of  Libby  Reservoir  from  August, 1 983 
through  July,  1984  based  on  30  m vertical 
tov\T  samples. 
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PERCENT  COMPOSITION  PERCENT  COMPOSITION 


Figure  20.  Length  frequency  distributions  (by  0.5  ran 

length  classes)  of  Daphnia  spp.  captured  in 
30  m vertical  tows  in  three  areas  of  Libby 
Reservoir  from  August^  19^3  through  July.;,  1984, 
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Table  17.  Estimated  densities  (number  per  m ) and  biomass  (grams 

per  m^)  of  Diptera  and  Oligochaeta  within  the  reservoir 
bed  of  Libby  Reservoir  in  frequently  dewatered,  occa- 
sionally dewatered,  and  permanently  wetted  substrates 
(by  area)  during  October  1983,  March  1984,  and  June 
1984. 


Frequently 

dewatered 


Tenmile 

Rexford 

Nunt>er/m^ 

Diptera 

October,  1983 

40 

125 

March,  1984 

86 

180 

June.  1984 

14 

4 

Average 

47 

103 

Oligochaeta 

October,  1983 

4 

4 

March,  1984 

4 

0 

June.  1984 

0 

0 

Average 

3 

1 

Grams/m-^ 

Diptera 

October,  1983 

0.153 

0.328 

March,  1984 

0.260 

0.756 

June.  1984 

0.026 

0.036 

Average 

0.146 

0.373 

Oligochaeta 

October,  1983 

0.001 

0.001 

March,  1984 

0.007 

0 

June.  1984 

0 

0 

Average 

0.003 

<0.001 

Occassionally  Permanently 

dewatered  wetted 

Tenmile  Rexford  Tenmile  Rexford 


126 

667 

312 

2,867 

256 

351 

334 

1,329 

319 

527 

628 

1.417 

234 

515 

425 

1,871 

7 

183 

50 

301 

110 

149 

52 

272 

305 

455 

158 

689 

141 

262 

87 

421 

0.328 

2.443 

2.252 

10.901 

0.820 

1.231 

1.485 

3.977 

1.261 

2.277 

3.555 

8.336 

0.803 

1.984 

2.431 

7.738 

0.007 

0.177 

0.025 

0.471 

0.203 

0.175 

0.056 

0.261 

0.224 

0.347 

0.203 

0.928 

0.145 

0.233 

0.095 

0.553 
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Table  18.  Comparison  of  total  catches  of  macroinvertebrates  from  the  surface  of 
Libby  Reservoir  (number  per  hectare)  in  limnetic  versus  near-shore  zones 
during  August,  1983  through  July,  1984. 
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difference  (N.S.)  at  level  of  prc±>ability  indicated  (p>value) . 
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seasonal  conparisons  yielded  useful  information  (Table  19) . The 
aquatic  component,  primarily  dipterans,  were  most  numerous  in  the 
spring  catch  in  both  the  Tenmile  and  Rexford  areas,  but  the 
abundance  of  this  aquatic  component  dropped  during  the  summer. 
Surprisingly,  winter  densities  in  the  Rexford  area  were  as  high  as 
fall  densities,  but  this  resulted  because  March  catches  were 
relatively  high  while  nothing  was  found  in  January  and  February 
samples.  Biomasses  of  terrestrials  were  generally  higher  than 
aquatics  during  all  seasons  (Appendix  K) . In  general,  the  orders 
Hymenoptera  and  Homoptera  were  the  most  abundant  terrestrials 
(Aj^ndix  K)  . 


RELAnONSHIFS  BEINEEN  RESERVOIR  OPERATION  AM)  RESERVOIR  HABITAT 

We  will  be  using  our  computer  map  data  base  and  the  GEOSCAN 
program  to  estimate  available  habitat  for  fish  and  fish  food 
organisms  based  on  reservoir  elevation.  This  program  will  allow  us 
to  estimate  volumes  of  water  by  depth  contour  (and  temperature 
contour) , amount  of  near-shore  habitat,  etc. 


IMPACTS  OP  RESERVOIR  OPERATION  ON  THE  FISHERY 

We  initially  envisioned  our  study  would  focus  primarily  on  the 
maximum  annual  drawdown,  however,  v^ile  that  aspect  of  operation  is 
important  there  may  be  other  corrponents  of  reservoir  operation 
which  also  affect  reservoir  fisheries.  The  timing  of  refill  and 
drawdown,  the  length  of  time  the  reservoir  remains  at  full  pool, 
the  problems  associated  with  an  inability  to  refill,  and  the 
volumes  of  water  discharged  are  all  operational  components  which 
may  influence  the  fishery.  We  began  this  study  as  a two  phase 
project.  Phase  I (from  1 May  1983  through  30  September  1986)  was 
proposed  for  collection  and  interpretation  of  data  and  the  develop- 
ment of  a data  base  to  evaluate  potential  operational  regimes  on 
the  reservoir  fishery;  and  Phase  II  (from  1 October  1986  through 
30  September  1987)  was  proposed  as  an  effort  to  cooperatively 
evaluate  all  study  efforts  on  irrpacts  of  reservoir  operation  for 
establishing  an  operation  plan  that  would  protect  or  enhance  reser- 
voir fisheries  while  meeting  other  hard  constraints  (MDFWP  1983) . 

In  developing  a sampling  regime  for  the  study  we  wanted  to 
ensure  that  our  sampling  effort  collected  all  necessary  data  needed 
to  model  or  predict  inpacts  of,  as  yet,  unknown  reservoir 
operational  regimes.  To  that  end  we  contacted  a number  of  experts 
in  the  limnological  and  fishery  modeling  field  including  the  USGS, 
Aquatic  Ecosystems  Analysts  (former  members  of  the  Fish  and 
Wildlife  Service's  Reservoir  Research  Group)  and  the  Fish  and 
Wildlife  Service's  Habitat  Evaluation  Procedures  Group. 

We  eventually  subcontracted  with  the  USGS  because:  they  have 
been  measuring  physical-chemistry  parameters  and  primary  production 
in  Libby  Reservoir  since  its  formation;  Woods  and  Falter  (1982) 
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T^ble  19.  Seasonal  average  densities  (number  per  hectare)  and 

biomass  (grams  per  hectare)  of  terrestrial  and  aquatic 
macroinvertebrates  captured  on  the  surface  of  three  areas 
of  Libby  Reservoir  during  August  1983  through  July,  1984. 


Number  per  Hectare 


Season/Year 

Tenmile 

Rexford 

Canada 

N.S. 

L. 

Avg. 

N.S. 

L. 

Avg. 

N.S. 

L. 

Avg. 

Summer/1983 

(14  tows) 

(10  tows) 

(8  tows) 

Terrestrial 

93 

463 

278 

23 

58 

41 

29 

38 

34 

Aquatic 

21 

19 

_z 

10 

9 

_4i 

_32 

Total 

110 

484 

297 

30 

68 

50 

46 

84 

66 

Fall/1983 

(16  tows) 

(14  tows) 

(10  tows) 

Terrestrial 

18 

22 

20 

19 

64 

42 

483 

43 

263 

Aquatic 

_ia 

JLd 

_10 

9 

19 

14 

100 

10 

55 

Total 

28 

32 

30 

28 

83 

56 

583 

53 

318 

Winter/1984 

(8  tows) 

(12  tows) 

(frozen) 

Terrestrial 

— 

— 

— 

47 

49 

48 

Aquatic 

— 

— 

— 

_3 

22 

_n 

Total 

— 

— 

— 

50 

71 

61 

Spring/1984 

(28  tows) 

(26  tows) 

(dewatered) 

Terrestrial 

58 

49 

53 

16 

11 

14 

Aquatic 

m 

158 

135 

258 

436 

Ml 

Total 

171 

207 

188 

274 

447 

361 

Spring/1984 

(18  tows) 

(12  tows) 

Terrestrial 

83 

52 

67 

143 

51 

97 

50 

43 

46 

Aquatic 

_21 

_2i 

25 

25 

21 

25 

JLZ 

21 

Total 

106 

78 

92 

168 

68 

118 

75 

60 

67 

Grans  per 

hectare 

Sunmer/1983 

(14  tows) 

(10  tows) 

(8  tows) 

Terrestrial 

4.37 

1.91  3.14 

0.48 

2.93 

1.71 

0.72 

1.41  1.07 

Aquatic 

0.12  0.09 

0.01 

0.01 

0.01 

1,02 

1.04  1.03 

Total 

4.42 

2.03  3.23 

0.49 

2.94 

1.72 

1.74 

2.45  2.10 

Fall/1983 

(16  tows) 

(14  tows) 

(10  tows) 

Terrestrial 

0.06 

0.20  0.13 

0.15 

0.90 

0.53 

1.26 

2.51  1.88 

Aquatic 

0.05 

0.01  0.03 

0,31 

0.05  .18 

Total 

0.11 

0.21  0.16 

0.15 

0.90 

0.53 

1.57 

2.56  2.06 

Winter/1984 

(8  tows) 

(12  tows) 

(frozen) 

Terrestrial 

— 

— 

0.12 

0.14 

0.13 

Aquatic 

0.01 

0.03 

0.02 

Total 

0.13 

0.17 

0.15 

Spring/1984 

(28  tows) 

(12  tows) 

(dewatered) 

Terrestrial 

1.33 

0.23  0.78 

0.40 

0.08 

0.24 

Aquatic 

0.47 

0.60  0.54 

0.71 

1.48 

1.10 

Total 

1.80 

0.83  1.32 

1.11 

1.56 

1.34 

Summer/1984 

(18  tows) 

(12  tows) 

(12  tows) 

Terrestrial 

1.56 

0.44  1.00 

1.96 

0.70 

1.33 

0.69 

0.35  0.52 

Aquatic 

iLia 

(USL  ILM 

(L2i 

0.05 

1.02 

0.07  0.55 

Total 

1.74 

0.54  1.14 

2.22 

0.75 

1.48 

1.71 

0.42  1.07 
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developed  an  initial  primary  production  model  using  Libby  Reservoir 
information;  and  they  maintain  a District  Office  in  Helena. 

The  uses  had  no  one  available  with  expertise  in  both  modeling  and 
fish  population  dynamics,  consequently,  they  were  able  only  to  look 
at  simple  linear  relationships  between  reservoir  drawdown,  fish 
growth,  and  fish  abundance  ( impend ix  L).  They  found  no  significant 
correlations  which  was  not  surprising  since  the  majority  of  past 
research  has  suggested  multivariate  control  of  fish  abundance  and 
growth.  They  developed  a flow  chart  outlining  a modeling  strategy 
with  our  assistance  (/^pendix  L) . Gene  Ploskey  of  Aquatic 
Ecosystem  Analysts  suggested  we  wait  until  the  end  of  our  data 
collection  period,  allowing  the  data  to  direct  our  modeling  effort. 
He  assured  us  that  our  sampling  program  was  rigorous  enough  to 
cover  the  development  of  a trophic  dynamics  model,  but  cautioned 
that  we  may  want  to  develop  a conceptual  model  because  our  sampling 
period  may  not  be  long  enough  to  develop  a realistic  trophic 
dynamics  model  (Appendix  M) . llie  following  is  a discussion  of  what 
variables  and  relationships  we  will  be  evaluating  to  predict  the 
effects  of  reservoir  operation  on  fish  populations  in  Libby 
Reservoir. 


Fish  Food 


Zooplankton 

Zooplankton  production  is  dependent  upon  primary  production 
which  is  dependent  upon  sunlight  energy  and  nutrient  availability. 
Woods  and  Falter  (1982)  developed  a multiple  regression  to  predict 
areal  primary  productivity  in  Libby  Reservoir.  Three  of  the  six 
variables  selected  in  their  multiple  regression  (euphotic 
extinction  coefficient,  stability  of  primary  thermocline,  and 
flushing  rate)  are  controlled  or  influenced  by  reservoir  operation. 

The  amount  of  useable  habitat  available  to  zooplankton  also 
may  limit  zooplankton  production.  We  will  be  documenting  changes 
in  zooplankton  standing  crops  related  to  changes  in  reservoir 
volume  within  the  upper  30  m of  the  water  column.  Thermal  regimes 
and  hydrodynamics  within  the  reservoir  also  affect  zooplankton 
production.  We  discuss  thermal  model  predictions  later  under 
"Modeling".  Flushing  rates  also  have  been  shown  to  contribute 
significantly  to  zooplankton  mortality  rates  (Ploskey  1983)  and 
there  is  evidence  for  that  happening  in  Libby  Reservoir  (personal 
communication,  John  Irving,  Argonne  National  Laboratory,  Argonne, 
Illinois) . 

Lastly,  predation  upon  zooplankton  populations  by  fish  in 
Libby  Reservoir  is  probably  the  largest  cause  of  zooplankton 
mortality.  We  need  to  investigate  the  effects  of  concentrating 
fish  and  their  zooplankton  prey  by  reducing  pool  volumes  on  both 
fish  and  zooplankton  populations. 
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Benthos 


It  was  obvious  from  our  data  that  reservoir  operation 
negatively  impacts  benthic  standing  crops.  We  need  to  determine 
how  fast  recolonization  occurs  and  if  recolonized  areas  ever  fully 
recover.  Preliminary  calculations  indicated  fish  production  could 
be  increased  by  as  much  as  15  kg  for  every  hectare  of  reservoir  bed 
that  changed  from  a frequently  dewatered  zone  to  an  occassionally 
dewatered  zone.  This  increase  was  based  on  the  assunptions  that 
all  Diptera  are  eaten  by  fish,  either  in  the  larval  or  adult  stage, 
and  assumes  a food  conversion  rate  of  75  percent  (Leidy  and  Jenkins 
1979) . We  need  to  estimate  generation  times  for  Diptera  to  deter- 
mine the  production  of  Diptera  within  the  reservoir.  The  relation- 
ships between  temperature,  depth,  substrate  type,  and  length  of  a 
single  life  cycle  for  Diptera  will  be  used  to  make  these  estimates. 
Ihe  above  values  will  be  obtained  from  the  literature. 


Surface  Blacroinvertebrates 

While  surface  macroinvertebrates  have  been  shown  to  be  an 
important  seasonal  food  item  (McMullin  1979) , their  distribution 
was  temporally  patchy  on  the  surface.  Numbers  of  macroinverte- 
brates on  the  water's  surface  were  often  higher  in  the  limnetic 
zone  when  corrpared  to  the  near-shore  zone.  Ihis  result  may  be 
caused  by  the  use  of  surface  macroinvertebrates  as  food  by  trout 
which  were  found  at  a higher  concentrations  in  near-shore  zones. 

The  direct  relationship  between  surface  macroinvertebrates  and 
surface  area  of  the  reservoir  is  obvious.  More  surface  area 
provides  a higher  probability  an  insect  will  land  on  the 
reservoir's  surface.  Ihe  aquatic  adult  component  of  the  surface 
macroinvertebrate  community  during  the  spring  probably  provides  an 
important  food  source  for  adult  Salmo  S£^.  prior  to  their  spawning 
migration. 

Pish 

Distribution  of  Salmo  spp.,  bull  trout  and  kokanee  appeared  to 
be  regulated,  in  part,  by  the  thermal  regime  within  the  reservoir. 
Changes  in  operation  which  alter  this  thermal  regime  may  impact 
fish  production  by  regulating  fish  distribution.  We  have  shown 
that  the  relative  abundance  of  Salmo  spp.  were  generally  higher  in 
near-shore  zones  conpared  to  limnetic  zones.  Quantity  of  near- 
shore habitat  may  change  differently  than  limnetic  habitat  with 
changing  reservoir  elevations.  Salmo  spp.  were  much  more  abundant 
in  the  upper  geographic  area  (Canada) , during  the  fall  based  on 
purse  seine  catches  (Table  8) . Early  drafting  of  the  reservoir 
would  limit  the  amount  of  this  upper-reservoir  habitat. 

Lower  water  elevations  provide  less  water  volume  and  concen- 
trates fish.  Hiis  concentrating  effect  during  the  spring  may  make 
gamefish  more  vulnerable  to  anglers.  Tiere  has  historically  been  a 
popular  early  spring  fishery  for  Salmo  spp..  lliis  aspect  of  reser- 
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voir  operation  may  be  desirable  depending  upon  angler  harvest 
throughout  the  year. 

Water  temperature  affects  fish  growth  rates.  Reservoir 
operation  influences  thermal  regimes  within  the  reservoir.  Fish 
growth  is  also  regulated  by  food  production  and  availability  which 
reservoir  operation  has  been  shown  to  influence.  Analyses  of 
otolith  samples  will  provide  estimates  of  seasonal  growth,  allowing 
us  to  compare  seasonal  growth  between  years  of  varying  operation. 

It  may  be  possible  to  evaluate  operational  effects  on  fish 
abundance  by  comparing  annual  trend  and  population  estimate  data  to 
past  reservoir  operation.  However,  deep  drawdowns  and  the  initial 
"trophic  surge"  common  in  new  reservoirs  occurred  simultaneously 
making  it  difficult  to  determine  causative  agents  in  the  response 
of  fish  populations.  Using  fish  growth  and  abundance  information, 
we  should  be  able  to  derive  an  index  of  fish  production  on  an 
annual  basis. 

Estimating  total  and  angler-harvest  mortality  rates  for  Salmo 
spp.  between  the  spring  and  fall  of  successive  years  will  enable  us 
to  estimate  natural  mortality  rates.  We  assume  natural  mortality 
is  influenced  by  reservoir  operation.  We  also  have  estimates  of 
spawning  mortality,  so  we  should  be  able  to  tie  natural  mortality 
directly  to  reservoir  residence. 


Modeling 

The  goal  of  this  study  is  to  provide  decision  makers  with 
information  describing  impacts  of  various  reservoir  operational 
regimes  on  the  reservoir  fishery.  To  that  end  we  will  be 
developing  a predictive  model  which  will  allow  decision  makers  to 
understand  trade-offs  and  opportunities  provided  by  different 
reservoir  operation  regimes.  We  presently  feel  that  a conceptual 
model  will  be  a logical  starting  place,  based  on  recommendations 
from  Gene  Ploskey  (Aquatic  Ecosystems  Analysts,  Fayetteville, 
Arkansas).  After  analyses  of  our  data  this  spring,  we  will  develop 
a compartmentalized  conceptual  model.  If  we  find  that  our  data 
suggests  development  of  a trophic  dynamics  model  is  appropriate  we 
will  begin  constructing  a predictive  model.  We  feel  confident  that 
some  type  of  predictive  model  can  be  developed,  but  we  may  need  to 
look  at  relative,  rather  than  actual,  values. 

Natural  variation  within  fish  and  fish  food  components  will 
probably  be  high  enough  to  limit  confidence  in  values  predicted  by 
any  model.  Regardless  of  the  type  of  model  developed,  the  modeling 
effort  will  increase  our  understanding  of  how  reservoir  operation 
impacts  a cold-water  fishery  and  that  information  will  be  valuable 
in  deciding  how  to  best  mesh  reservoir  operation  and  the 
maintenemce  of  a reservoir  fishery. 
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We  plan  to  "filter  out"  much  of  the  environmental  variation  by 
using  a thermal  predictive  model  (Adams  1974) . Since  teirperature 
regulates  or  controls  many  functions  of  both  fish  and  fish  food 
organisms,  we  plan  to  predict  the  therml  regimes  within  the 
reservoir  using  the  above  model.  Tliis  model  will  allow  us  to  hold 
environmental  variables  (volume  of  inflow,  teirperature  of  inflowing 
waters,  and  solar  radiation)  constant,  while  evaluating  iirpacts  of 
operational  variables  (discharge  volumes,  timing  of  discharge,  and 
depth  of  discharge)  on  the  thermal  regime  within  the  reservoir.  We 
can  then  use  these  predicted  thermal  regimes  to  predict  fish 
distribution  patterns,  primary  productivity  (using  Woods  and 
Falter 's  (1982)  model),  zooplankton  productivity  (as  an  index 
related  to  primary  productivity)  and  losses  (based  on  generation 
times  and  flushing  rates) , and  fish  growth. 

We  can  also  use  changes  in  reservoir  surface  area,  reservoir 
bed  area,  reservoir  volume  by  depth  (teirperature)  contour  interval, 
reservoir  volume  by  zone  (limnetic  versus  near-shore) , and 
shoreline  length  estimated  by  the  GBOSCAN  program  to  quantify 
changes  in  habitat  used  by  fish  and  fish  food  organisms.  Using  the 
GEJOSCAN  estimates  in  conjunction  with  the  thermal  prediction  model 
will  enable  us  to  predict  operational  effects. 
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RBGGMIENDmONS 


1.  Continue  sampling  as  modified  below: 

a)  collect  food  habits  data  through  August  1985,  then 
switch  to  either  cursory  or  pooled  food  habits  analyses; 

b)  investigate  desirability  of  sampling  macroinvertebrates 
on  the  reservoir's  surface  more  frequently  during  the 
period  i^ril  through  October; 

c)  use  seasonal  rather  than  monthly  horizontal  gill  netting 
to  sample  fish  abundance; 

d)  collect  three  zooplankton  samples  per  area  each  sample 
trip; 

e)  relate  euphotic  zone  depth  to  secchi  disk  depth; 

f)  sample  for  burbot  larvae  by  towing  an  icthyoplankton  net 
during  March  through  June; 

g)  electrofish  representative  tributary  stream  reaches  to 
provide  data  for  developing  a model  to  predict  potential 
of  tributaries  for  producing  trout; 

h)  change  seasonal  criteria  to  represent  thermal  changes 
within  Libby  Reservoir;  and 

i)  conduct  diurnal  sampling  of  both  fish  and  fish  food 
organism  distributions  in  at  least  one  area  of  Libby 
Reservoir . 

2.  Conduct  a creel  census  and  economic  survey  of  the  recreational 
fishery. 

3.  Use  log  (x+1)  and  log  (x)  transformations  on  all  fish  and  fish 
food  data  for  statistically  comparing  annual  abundance. 

4.  Sample  downstream  losses  of  zooplankton  and  fish  below  Libby 
Dam. 

5.  Delay  trophic  dynamic  model  development  until  enough  data  has 
been  collected.  Use  present  information  to  develop  a 
conceptual  model. 
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